


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1984 


Normal approximation for response time ina 
processor-shared computer system model. 


Pornsuriya, Suriya 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/19471 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


f (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

| | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


A hee ee ee te 


? 


=" s9 


A 


* 


St i atl 


“For = 


| ee 


a 
— 












































‘ea 
2 eee il oe el A Patel bee (hr 
Se Sey ee ee ee, he ees Wha ON) 
Ae ee ener er r wa eS ley TY Lt eh Yh eae AUN ROT eh bs tet va 
.f et | a et one Le ee Pn eer Lh, A ate et Aas ken Yk He Slice’ is 
[nF a a be by ee Lay thee te en | Loe ee ai ot nee WI Ve teak pe Was Re OR Bi eteeat tk Coe YY 
. i Ly Reg) Pah Ae ML a o@a ar i ar Safes tas 4 5 Ep Srey s ee ee te | bed Ae ae ard bhi le kw tS 4 7 7 rete TY 
5 eds Cle ye) | pie CS Ty oer tae ‘ Narra Deed UU AG ae by eh ee et Ae ol ee A pte de eth ay 
u Ld PeLen rn oe an JI Ty fers te ae L oe he oat Bate Pia: alah COL oN eee ee au Chea en ek eo he " Sek hed hae oe a 
iu oO ee pale rsh y Lu Jl nL ar, “ae eer er i | oem og Mb y —" yt 4°14 8 Ba ty AF, Mik 8 any! a Te ee L ec od ee 
wo. 4 on zl Jetty iene met SNe iy a - aNd " 488 Mile ee ay a, te asa Ke oe ey itty added s 
1 pi be ee Ve gry Ps o . eC ate hu Ue VA MY Le ae a eet Lo 4} Beta: ee ee (cienkon a ee 
1 D e . aw co ee a ea Lee Oc ee er es se ‘hos A Ys Kd ave en te ey ee ok ok | be-Tiep Mt tele. oT. Ty 
q ay eho ee tn, aa rte ona ae Vat ges) Cra a mes URE ae Pe ee ee A ce eh eS fn Peete thant ee he 
ry Ca aT Le ot em Ered as o Pan] ‘2 aan ee CS i oe at) 4 ra ee ee VR 48 we Oe ted ; ed 
ql ‘ . t “ray a4 ee LP ar Sera ae SOC 2 c 89 ae tae ee ee Ph oe ee Pre Tee et] Ci) 4m Ob Ue gran rogliy Sept tanh 
ry ‘ be i me | | Sart nen oq Deo) OGRE Va ue rete tag beh Pech He ee Sey a ee oe ee or ey eee ee « wie I tata ch dean be lead 
Py F 1 oo. "eaten ’ ee . + saat ry e a % ory cy ey eR lela Peet So Sane Cara 8 ew ae eT hat POR eK ORES Se ‘ Lok 
a er bh ae ad ea Uae. 2D oO ia Urey Nremwa owt, ws OVA ey berry bare teh Me Bel eee a ey 4 ie Apt ttel teal ty es te tetheriact 
ra Lh OL Te Ce Wee rs een Sara at ear yet Wie et Ce PL ee (eV ke na Rg ee in aint Yay LAS ate ree 3 re OO ee md bate Keates Adal, 
O ta 4 u U i ee pe EOE ect ass Ce a er ar Le oe Far era hay tay. ary e's “24 Loe aE. eee SNe ica ance rte egy tae ey 
a ee ee Ee LW aU eee esl or Viren a ae at ary ee Ree tte Ley Ser to Le Li Ed Maeda AK, ATOR er wt A Serlghindgyst, oot WU ae Ree Mette th tis 
q Ca) bt INE me re Hur aa | Xe ee ee ty ta MO ba el ST te eS er AWE o-9r8.9 4 beled ea a Gh eiehirde et NAc test tt ti hee 
? 7 Ly D -* LO, eee fet Ye PA YL Uae or Le tee ee Use tan mes Pec Pd De at 5 AAA Be, rig hr keg Le . i PibA Vay bids Shot oa PLY bhatt e  e 
f . Oa C saw ay UO Ee aor ee UE eae hy Pati A URN Ea irda PAM AeA var beg & er de wey con ae a 4 i! We de] Saas eerie | Maral Nerkseetotee date eee 
te fy ry d aay | a La iN es, she Mar 5 eee er eer oa | Lie 4 Sioa ie ee Te eer Rea oe Oe OP br Se o >" tere ae AM Mierbehat BF et bee we A: at + vs mrhabet Ue pecdicn 
ey a) "4 no er eet Crs Ca ". "qi i. a Oe ok ay ate art A ce i te ad ale se AL hd Dt ABih cL re phe at re ipl Note Ee ee 
o Pan LoL eee ara i Sra cr ot ery a ee rr ae rh We Utsrn ley kk MLE Oe YE te hr ead A ee Bid ek) aa erie aa ed sy ae Re pha aoe 
so 7 a) Or ier Cae oT Hi ‘A sa, Ae OE YR er par Pale a Che 93Woe St ae ay EP shy s Ayrye mame Ake dzanyn beh Wh) Dolebha hOk  s tn ded 
) G : ' a oS TLE en a LM ar “bUe Med ae: Aa eM re hal 39 ad Mme es tLe beat he baa a Cr bie GE Ok Warn 
a] * : a, a fea Cy a | “ te . Car a or a a) Lie or ry ' an a oe ay te | ay be ets ao ew hs 
BE o a a t's "Nya hi Sd ate , ee “oN Retest ep & mz baa ee Ses 
* tes ee ts 1 | i ee se ae | a ee ae ee ee DS es hy ic lM Si areas 
ae A Testa aus Yor) i SoC A | btuesy cL, ee a 8 OO en Or Sst len es ta 
Cr] ee iD L o 1 ar] Rory ‘ (a bwey 4, ATs 1 tg a Ia oe a Baselea del. oy & 6 Deities ee 
t u ied ie | Ca "~ebuen Vee se td nog vo ota. batik ie AC cr am ale kL are arty up es olay eli derteth ate 
PT 5 4 L elvan A : i he ALT Ee aC 1) "XN age ' ee a Pe en} pee Pute de lth KN oe 1 Aas 4 tS RAD Mey wh ly hither ok re a Toe el wr gets. 0 
' nig ' 1) 1 4 4 t Da ca i at é Pry bet pam tee ny re ee Me Urn ey Sastips lea VOvab ay obey LOR ee eer me 4 ba 3 te ten bathe re nt ea ee er CA 
ae foe "a o : 5 L pe ie Reo nner arte or Lo ee Tec inariren Le Se i a Pe ea Leh eC ee ee EL) pe Melia det dS or tim Nor A aa tet ees bw 
sd o "4 as re ie EC) fer 4 *, Wave vy yn bp a a Wa eee nS MAE As g he ef ees Sy i Arotel th fete oy Shekel Wel Sb Pied feotin Nhat Send a 
aa u a aa 4 ® aL Me GETS, Feta pte DR Lyre Fac Vor ay dot nah egy AMPA LN on ey Sieh CN hey or ater) bg bndnan eee ty ellis o 
o a) ve 4 4 DC ea eet bs is a ei) Seman sy Oe ek eT ae 4 ae | ya bata Sh dl Cee rere Ay iy Am? nt tee Aart Bement ten Cle 
A 5 u ‘ er "Oe eer eer a se ae Oe ora Male igi bl Se ire War) MMs 9 a ae Aly Le es ay rh bt of acaithalra a htehe ene 
"4 an er ae ‘ veye oar) Laval har ser rat a] rn aegy ary tag Brluy | oe ce ay ta Aaa ae BE Pir ve ir ir aA atne CT 1) a ea oy wah cael ee rey le vs 
‘ a) ‘ aL) Lea ‘ r a a] bene, Le ay er Te 54 5 MA he, Ssh RY my ee Oe ee See Ee a oy Atay To aa hone bt Tn ha b Lia tte ee es te C 
7 5 ee ee 1 <P Lape an) ‘ $8 wy sy cee Kees UN ea Oy wee anae A ea, rc YY A A ire Ld TOT eet tits Mibh Yeo te re 
ae ert a at) er awe ys earn nae 0 De Et tip Smt VET Er ent UE RE is bd MU OO Or Hr Shar Hea RIAD bay Sahni tae ot AC ne beh cs be eA he ota hh we Ainter tech hhc 
’ cn ry an) a : ™ Lae aL at he ToC) Pa on a ar] n ret Wr Wer Pes i Vise enim ae 9 *“Wers Te ae eke Lear wean HQ tum bho, A ; piel tot) behead ik Se ie Tne ae oe ‘% 
ary i mn Rohe re sae er et NWS) ey Aa WeeeWpar y a ety ato vat i On a whe en O50 ca eee Er bial et 4 bit Withee ty ow, hi-Aaceh: cow, ae 
Cee nl Ooty Ya! "yon o rere Terie A ho Se bi Te | Vetus we, oak edn La BO bot i Sa 42 WN ng hah be Mei dy td ce tit eto ee nin 
g i oo rg ‘ Le in Oe Ty Seems US Toe ah ee Br ery a Re CW er aha ra OY ares 4 rr. rh ae bgt Mote A Malan ohn Darter er RA earh hs 
a o% 1 es poe What at eer ' a Per Ln) als Te pa 4 n aaa A Se ha my ar] aor we + ee eee Le he eo a wey hk a PMO: Parka Ate, leat te ee ee 
= ey i r ewe leas Taree D & nN Pn er ar) ary Seen ® A ary a Ph LL J DL ts a Lied Ok We Se ne aL a) etek is ee er epee oat ST ears 
pl ae u ' i en ee " a +a en 7 Lr en ee a ee Very mae gua Be eC CRE RC ae ie Wa) Reo aera eC Caer eet 
’ f La : er ey ar i at | rans BY rece rn ova ara Y b ae ey BAe wren ow Lh eee Oe Pe ice a) Yet ended eet ey helinde bona fe ee’ 
O A ‘ 7 O Lr ot ry Un OGeei starr reer er | Wyte n na + ry ee ere Pr Nie oe er oe rd #0 4 aoe or 6 ee tO ee ore tchae tes be dodiere 
oO ) he i | hg uu v u' 6 eye aCiey ye rh ie MeEA erp ee any es c a Ae aL ae a eae hee RTE heed etc. UE. PON ae What oe ene 
Ca rar o n Cat ame tue Cory ry Ce Ae C1 SERVE Wks, Pera gy 2 ML ea i Ue ans PT Cie eee, a Ur ru eect ee ey Stet teagan atod, Taree 
' 1 Py "4 + ry Pare) ’ Un ery Pr a} [ers ha oh fom er ee o ie 4 as i ha a a rey oe ee ae et be aeb) ute. Wa betel ee ex Ly Me tert tok ee 
q Cr i er o OV wi van a eRe ELIAS). Sere rye re ar LOST a rie ch year LA am At eas Mia hk AY A A Gg FL A tetting cept Ge Ba Ec A A LT SE lode ee 
. car) Cn ro] ar] my 6 o4 oan ry brea ge oa ra oN 18 ae see a Le ee een aero “%%bheengtrve bch eee a Cabal UEC WT ae neat Mie Mne th iat PS 
‘ae oO ‘ eee “ 4 cf "wha polis Lae ST ry eta WA cs ae ‘ t. 1 ee arses A ah REO ee Ca, ba | Me Abdel sid Is ghey ty i De ee ne eee me wy bubdos ke oon ee 
a) oh ay cary 1 Yi We Ye eh Yat Hr a ry CLO! Sarina mw 4, re ee Ser as Oe as MOM Rythme gt ke be ae 9 Le a ya bd ak) a Lk ie ouy % Wari ete Adare Ke wo 
’ ' ‘ ’ ’ . ’ on ce a Ta Cry . ‘ on wa os Pers OY Moria. % tr body ce aL Es J oe o 1 a i ie ee LIT) htt at el darted ee ae Se 
Laat Year | bY ‘ ’ ’ . ny ’ R46 cn ey Oa Pht rn ' Ble meas) eC ee i an cee Cs ay bial et ie Pe Wan fmt 4%, f Was er betsy id ey Ste Sale eet Oe 
' “so oe 4 Cy pew ee trier ny at NCAR Se ey ale Os ala ayer) ry ui Mw abag ie, ce ae oe SL Sar ee ie <» PARR ae Thy Metal) RE NI A ity ets 
a4 iu m4 | bs one ee + Q " Une Bee Se YO Sa ese Sar Har ee. SO oe oe ee ey sgl > eM Bah Je ee ee UR ack es ters Lt ee nt oe rie tuto eee ee 
i ' hoe * © te ee ee ot s : | Ln Ci ac A Oy Wey ay rn ar were FR Ny one § \ ht eee ny a Ie we ha Wiha ack oot Tr en ta 
Fy ary . eo be c Pa ' ar] ee bola Wace eas te erry pC ae act sey ane v ery oo an en "te ee or yon At Go oh Warant s “Rilke m, Ue Sa COR beri tete ties 
1 Let aaa a) ‘ a i een hy oO an Sea ee Ser ces err Shee eee Ay ne + SP eonye we ba aot | Bh Woke te her ee ore we 
r ’ r oy La . ee ri “a1 Liana? in boar Soar Corian Oa tal Or Met © me og. Hehode, Tagen see cre} “@ he 7 he eee aon AG Ee ee ‘Ter Sen gat 
te A Fl oO bY »ad Can D i Ta Yat i ar) rt Ve oom 1 ry Lee Ce. es, ier ta etl tong eR I ce area re dale kOe ere se 
. Ory ’ ve b CY RHea, aie Cr hit een aren con ee ers rary er ae ny a Pes ye SOX btn bing ts Las a te tod re) ith be be ny eee 
. A Del allo hey ga es ay - see ee ST or yen pi Saber d var {weet Cy iace v4 LA aOR err pat a, fe ta We Weer ha YV PAWN VE Ps ae oe pee a geet te Loner teres 
ae ar] ’ A A Fl a ' eo. % ry ry I EThon SA rt a rye EF ey ne) re aay yewa’s thats, 2. Se eee be beet Le tine} Dy Sele te 
rl a A ig Caries ’ a A rey Tr} Sar ter CAT een a ¥'S ry RORY Gales rt 524, AY. Os a A ee ee Mea ye \! Ce eS einai a Datel th tee 
ae ame 4 Vas Tiny f as 4 ie Sn yar cy ‘ SE ern oar, 2 aoe ee ee a ae ee 5 NY bel taal ra hr ORR Sch si eee 
0 A Caan Cn) a iL Wis an ‘o SNS ON Gea Lee rl a wi > Oe re 940 1 Men enn S ae en hey Lor ark a My a oe asa e rs pede Ee 
A A eT A a me oe er oe er oe RRS VG vel ir var te ty, rere th Siem: IF ict iy uecrs i Pa tel a a a a Pt eee Petite tanita Tey CotT esa aoe ee 
O A i bs a er Yaa cary | | a) uy ir ay {bee ye i er een, a ray er ir Shea sa Be Te Lite”) ae! tao dh | Dah eh ee es ee itt a 
Va ‘ D U Liak o * ae $ bene Ha her in : “° ci (8 hee ge Varo nt PO no Ce ea re es thee A a bee taht ee 
’ ry Cn a | ' a4 1s ’ a: ry ae ee UC eae a a I Parry 7 4 CU Pan aan Lan ph a oe ie aT, inte 4 hy ea a ere OMe re AR he 
1 Ld U uy i : aa ’ + ® ‘ ry ‘ ‘ . ag ‘ fa Lee | (4 Pe ee oe ete ry Se Wee ah an a a Aah ons Gir ssn wee We os med Lenk Ta 
oar] Car) a | -& Nara ee WU lent Vary ary oo ad aa aries ae hk cnr ra Pe 9 lead es Bete tot hy tang * Vana bene a ee belt fe 2 te ee 
i Lr ame er 1 eee) f Peay BN oa bees o1y aera chara » 48 1g SMV ay a bere ht ae. ek wera eee b ee Stee aE Way bach teak open ae by ta Wa ellen 
is Pie Si A ONG See argh ats 1 ¥ s ". ag en er ANE: ST Oe ape aa asia YY Ue Gr Vice Sil Ph oo ee Ue Te a tet ah betel ee Deh ots a belch toate here: 
oar) ’ ar La lr Cr) 4 "et «ff Cerrina ms a 9 can erty a a eres > ee S| MB Ald ee men ESE ble anday yt pan wr hoe bibs eo eee oe 
A aes ee ae tay fe nar 1? RUN CIioie man i ee ORCS vce tebe Ve te ate A ee ser ere et bee Ley Posie RW See NEN eee, aon trek 
“ D 5 i i Ce aa ae D Uae Pan a Loti aes marin ws Pa ee Sry eis veil s° der Sore RO MT iit) oy Mas dt Ta ee Se Ce eee 
4 o ) a end . ' Y gee 0 ete af] el ee 1 NS tas i ia Oa at We ORY “ Ne tel a Vaden Trt) ee a ARG Lat er ke era eerste Merb: teak 
o : E * LUCE VAL CSC mens ech Rarity) pies Carreras NO arg ay ea ie Sree) ae A Ie Ok bad a, ee To bake SOR ee ee eo eer EN 
"a Cary ary iD iD o.8 + eeu f . o. Cr ar “rea, CO arity be Tai. TNs gle ae Ae aS ST ee pins he. Goo IS Mn hed ahr 3 Shad 
Py fs fo a] . i . co . ’ .Y ns tr eae ma Cr i et Wa lage ty RN a) AS eT ate 7 ee ee NM Rae Cor eT aa) 
‘ iu Cn th ” ary NNe agar a, ee, Teeter od WR ist a wu bet ran LoL TOY oremar epy 0 cat rear ye ers JO a aie a a ad Ltr) Hon hate ta elk ae th] L ae Leh ee ee re bade Wh a woh 
n D 1 fie AGI hs o 1 ¥ oa oO Ys er . Le ae inte tey aia a bald ot ody oh ke os th eee tet SS Steffen. © y bbl Se 
e A ry A A "oe Cr | ay a arn pie Wy CY ar SA aa een Gt a Cary Phe 4 Ae dane een We ee Leh ae wee 
. fy A er a ’ "i Cn ar] a ' U *, " BY ry hy Cue eric vers '"aite ma. Sea ke iittadirdr Ak aces ct an 
oe 4 n a} me TL ee) ie toe . rT feu List * o bday “8 ony Sa ea te ae ee ‘eer i) 
‘i A i. N Poe oe (ene seer Tether Fae ener * He AED os F Sere Wear es Siva: tri Oe eee h 
‘ A i o Mase he ’ o . 1 va - : ar | toys ar ar] oa . NWO ae beh er a a ee Le a Ca ts a tee tek hn 
a L n o pl > a aS “ee wah Ly Orne ra > Peet mh "5 tase! a hr Sth Read Poe Lag eaters F 
‘i , ! Fy] . ’ Oh Ce 1 ry A "¢@ ’ Le Paar a i ~ ers Le Au ra ~~ ar arr ae “yew wae nk a Seta chee) an ye ar 
’ Ln rerrecy r ' ’ Oe 1 oa ; a er) Pi a " » ares ow ce ry ei th Lot Py en bar ee Pn on i iar) ada aly Ps ie ee crn a 
Se Py ‘ an wae A A A rae ery D Ory , he ny Fi Se earn ares . cra a ey ft i hn yall wr | PTY Ste gi hee eae yp oo 
i aL Lee ei o a a age OD ‘wou ‘ MWe Nett se hay ww! ef alotay a fa we pn Tn Leo ies or ee & time he | a 
oO 4 PF 0 # ot y ae? ‘ iy ev ws Cn ee ee he sy 4 ’ at cafe 4s A, ba as vi oe bs a7 shee ee CLE fee tel i PAB 
ar i * 4 BU ioe ce ql my 7ty en ne or ar . D ae Oe i 1) ee ee i) i) MA to ee Tee ie 
Pa ' fn % b- ea J Py Cis era a a ore ‘ere a LU on o Gis oe agen te yh Pa tes EI) eA Se in a's, Lae 
A A p " D if A eee Oa oy Maik ed ee SUCRE Ye 5 0 aE ad a en TY a Pats lathe Ceca eet tril 
he ey a) i ue er eee) hee Bee 202 ely 3 gS TN hg ek ih Se 
bi 0 > rs A han) Hee st Ces met Uk, Sane ry es aan eee ee iat) Rate hin ee erie “ 
Ay car) oe a) " ‘ car) U mNeny pet ele eo atirs fi Be Soh MAT Ree So Pea todd iL bes ales be azee =, 
0D 4 1 U 4 U n a fo 1 an ha ary . re . % vn? Bean Ur, oe Ee er are eget aa 
eon i i ra e Ce ry re | ryt a Fiat 9 Cn ae Ser a te ra 4 me Ls 40 Se tg ig OEE wy ae baa hed 
TT 4 of ‘ wee iu : n ot Lia ea u. aD a hes ry ra, ae iz Cte errs a et oy | "4 Meek Ra ee are ee 
1 rY 1 ‘ P an " or] ‘ a ry a ary 2 ea br oO a eat!) Saecr a ie - 7 Ne AE Ag i ant hae a en we . 
. oe bared i D Car hay yes Sa bt pont) rye O in Ue x A . ao 4 Camry Ti, © ae aes A Ea 1 hs % » ee Rey Cyr Tere 
f D Oo me a bi 1 A f bY Ve ache H Sr eee, a | oe oyu * di nee rial Wa! va vot ha he Dee at bs ET cree 
ri i . a ’ Pi aa! . bee eC TE ac eC ot » i - ‘ey 2 Tr oy a’ bn a Pets Lio bil hE ee ee " 
Lt r ry A = "eu i bd * oa Sa] ' hy i) iy ry ri Sy ces are ' a) re oA 1 a ate and eA ae Me cee A 4 amt, a ey St a: Fy pat? 
» D Oo oO Laie bs i Lie Brace iy. ar Woe yas Certs a a » 4 ons b} se et ' AP fie: o \wrnd i: he ee = he" 
Oo 8 @ anne . fn Oo a n ‘3 tas o Lea Po hd Py near LY 4 Fee § aa ih ok Wot We ay tits | a ee Tg) en oY a ee By ee 
eer ry . . ae lee, LS 3 Ie Dit ry Ce ee ers oS ay weet hae iid Piety bor Se als feu, su tat Wo. ho, he Lr a ies yo 
ks ry Cr rs B Di a it ’ ‘ oo Os bd rota mg wer Char) bd TY act Fa ry &s we™ 3° “a9 Pare ty he ar ary Ty gts yw Ue CF abd mer 0 iy oe Oe ee eed 
re 5 + ‘ cde Lor) ar A » / + at a reas ror ae 2 ae BLS wae Roses Gite ees ft Ae, Me’ a ee 9 Ss fe A 
i a a} ‘ Ld A . + 4° 1 an] Yet ay 7. re ey 1 ae <a nek a ee a Oe Ly i ar boeken Te “ag! FoR, idle “eb Aa EN OEP 
. . + ] oD Oo * * . Cd . u aC . ary Fethe Ta eg 08 4* af rd ar} PL Ca J % t ats 6 i iar <0 Monee L ea 
Cicer rs ny ry ba) ’ te ‘ ar) an] A LS. ae Delete 4 a ee + LP as eo tee ad gh ty Ce ou a ber) bear oa) a ay re baat sg Fe ei ate a 2 Ce 
FT Cha ‘ ‘ J ‘ in : ae ay ar ‘ Cte PS. Ps Lh eee a 7 iM 7 tr ° nd eg me Poser 8 a ee | at oral} ew note a Stee Le a ee Kk “> 
‘ ’ a har) iu i 3 teh * . Oo eat Lr 1 x, ot ae _ CY ~ yw “Pate wy Yee oft, gen, wahee me FNS eer ye 
a a] » . ry a eT) A U . ys Care rey k rs cre ry a) ) Pi wn rag cay a <5 ete Ne ey tee * fle. eet 
O ha + A ce p oan ame Pere URAC ieee TAO ae Sc a a Le) ea rte ae eee AA eee Lay a Cr ee. he 
A ’ a ’ 1‘ ear ry ria ’ r ‘ . | hy fe ri os r fo 8 Me r. ' ea ed Saf 2 enki ale By _ a ee ae Pe ae 
£ t rn 4 u te eee’ f Ds. iat ene tee ice erie) LL ST Irpeh ees * * f A * D re "atm tiar. “sh AS Vcr re fs 6 Nat met Gite, 
‘ i Laer . o % »' . ae L uJ La + ’ bd 4 - Thin oPyage in ee ye atrte th sn re) Sts 
ns rT ne yi OD a ry & * Fy ie as: Rat re <0 oe ot oe Soin 4 Po Net py 
Li . oe + ri . ny f vee 1 bee as he ew ee ae Ce [ir ae + od 
o : 1 Ld oa f oO A ra . a dy Sere foh 6 aye a eat ed oe f ¥, : 
u A UF n Crs era ers Meera a Ce eet 44? ¢,. 
D 1 ee ’ aac . nk ae La ry a at LW a meee ry 
2 , ee i rh Ls he er eS: She tee a 
F ie ‘ ns i O Ora | oe, ee Se a ae Dee a ATT 
' i. = s ees wid F Aten Et tae 2 1 arr, 
r ed . + cite Ce er ey U bs a Oe od he i) an ers 
', Lamar tt aay it ES Ac Sh ae .F “i 
Sal i Ly + rT as oa Cae ce ef Ls i JO A ‘eae 
i o fi 3 ny e i ft oy Peay 2 a De eee See otra 
. 1 ] o* . ad rs a r pen os “ee 4 tay Ke AH Oe eee 
. ry * rrr = + ry ar] ¢ ry 5 ia Fons s My a) Pay - Ak ae err ae . ca at) , 
Z + » i hi : : Lt a sf . ED aw a a ae yi Ca - pn ie eee i — a eee +e 
i ' so : A gee i i « ye an a) G ee Lt ae ee isn Pais Ww ° fe Pres i “a al aoe ser is fate & 
, ‘4 D oa , 1 or rn an ° o cee CRP eee im oh ie SL oy ae 2 ee 
7 y aa A slot eo geo 3 ARNT ial oi AN my TEAL S38 4,, Sera A a a ore Ln v>~ 5 la 
a A 0 G cn oa a * a - Ale at Fa Met (ect en PE peeat a ees Oh Peg a 
‘ . 1 Lan ae Wee Sh ‘ Pe ee « #heory r ee ii: rai Slee ae Lhe LP ay Ts, og Teen... hee. nih Py okt oh q 
0D - Cian) Pr) ene ete ’ ’ + «© af @.! « re hy + ae Pr eas ie Pe es) Pes oar Ah wis Le ai kare 
p oo teed Me can ; PARLE Shee ee ct eS era Bs NY ne eo eee Bri eee 
‘ a A G : Ae ee ee; , : t vo ae netae Fi H A A I ee Pc nee 
iY cay r LU ‘ i plete 0 Oy tae Pita ie tees "Sta fF Or Le] + Phy | rea Cheers 4 
. Fy Pare A A 4 Che a) ’ ° Pry ee ee No be 8 8 oo, , Aa me Oe Ty a wane ces a 4° Ate Lame scr re a 
A ar] ‘ 1 a) . + P ry er mer) Pir rs 4 Wa mq ooh, ny ema Ra Tees a ee Con yr eh 
1 1 0 Ca i 7] ’ nl Cat os A oe rm f %gt @ yp ae Li a ir) by er) corre at Ce OP ae ee a “ee, Fi tee wae 
A A oe 7 ’ ry ’ an | 4 o\ ’ a r a.) a LAE er ey a od Cae ae Cae ay ee Laan) a Pre Pha oa Me Tale es my 
. & i] Sa cy . 1. ’ Ly O 4 Yann a *Le a had f . ar) ne ' ae ode] on a a] mete 4° te 3 ret 
‘ D , ra ae = + ‘ ys Per aces EON .? . al a Cd Orin ones 6 . a, i Ce 2 ee ere ay). 
P pl a A h A a er) o A $ UL ae MeL Wren AM en oe Cir a Or ae fr aaa Ce ay caer ae eel a 
o ‘ 3 7 o 7 ™ 6 ' i M Cer A Dae Cr icin al Nee 7 a + Log LAGE Ba PC eat y A ar) a es elt dae ae | coeiat a 
a as oy D U 1 ri re. Par ha ria ‘ et te reer Par i oe Tr ar] er . rhe Late ? Nr . "em a ety nid ak Phld be ed | 
Cl Ca ee ' ' + Cn 1 . ae rn rs + i orn ’ wo Brae on Le Ar re a ee eae rear he LG Can 2] SHAR Vets ae ee) 
O ; rd 5 as PN eg crn i a) 4 oe ¥ Pia Dey te ean Pcl. a a a aes Le ar] Careers oro ogee a en Or oe er) 
ae ad Dae : D ecient  ) ‘ MALE CRC ot a5 Bee A ete OC ara. RAN haa of, Pott Ae Oe hair yA gee Er 
t) iM . i ies sen Pa td UJ tog i rs CSL tee? Wersaes '.' & 7 rir) 3 1 mae’) a, gieetcal S eee ako sXe eee) eet at, 
rs = i a ay . ‘ 2; UJ LT Tie eel leat 7 6 . . ) + ee ee gw, can Cn a a ta Lr ey ae el at oP L Se ok i a A : : 
n ror o Fr 7 ri eo A Pre ’ soe) use 6 4",, Peery f- OT sigan ee ia Le Pas tr fF, ° Wat and ae 
’ + Fd 6 . ar) . UY ta er J oad er wer Pi roar pet shy (OR 4 tg ary Caan a Pare ie) ee oo; Ce ae a ae Oar + eee ad ae Tr 
Ay AD rn) iD og ‘ u = s+ 4 agree oy Ur OO aC ae CC Oe 2 Oe) th eats a oo SY Re Le Cie ae al ; oe 
rs r ra rs rary Pr a er) ror 1 Ore a ee TP ee Year ro re cr ie La ra Pry Pre Pay ok a at om Cea eR a etapa? Codd aed tr 
o + ‘ ary na ‘ ' 0 ek Led Vr 0 geared EUS oa ERD Sor) nj) acre CL hr a Ye ot & a ee) Pee aes © We dit eae AM ree ss ore 
oo . A oa] ble liar Tat le a PY Sacer a ee es Cree eee a a 2 see yer ier Cadet a alt 1 oka) ere f oy ae, - gre 
OS Sy Gal 0 a ry ° . . nr | Ca fae eer . are ra nN an ca) ii Me Cee er) Ge ad Risers ad Ce R. Aa pr be ye Ae Dodane ee re 
f ar) A A ar f 3 ‘ ‘ eae Aeron ; ri res A P+e, ay SECO are ey ey ae ote yet dng fuss. OC ae Pe rs rae Ce eel at ie - 
+ a mo any ‘ ee ry Cee ae, A er er on Cia fel Lr Drees ig pe er ne a dL aa A Ow oy r= ms ited dl as i oY a ey, 
. oa? hale etd ’ . oo. ’ | . . ones ae Ji ee) er a fies got ar] Cale ae et Or ae TP Cae Lay) Rat a ee) pee Pat lee sae 2 Saal ee ad 
Cee art O D mee deme er Yaar Ya | LRT + yar A or] ‘ ee + mt ce up of # ay ee) bie ew ee A ed nc ee SR See Mh ew 3 fine 
¢ oan Cr) an , ra UeeCuraery Ea Tes Ca Get bier Mente Sito. a LL GP aie a ca ar) bk acd Sa arr) CR a ae ae | htadlbahles dod eh aed os | 
1 n "4 Ln er ae ' oe 8 1 ar) D Co ees ee od Cae Pe siar pat 1% Fy i lea ae rp lee weds | 6 wood Ye . dt er ar} ae Oe a 7 he ate | cpl at 1s tae Cay 
Ls u ° i ‘¢f nas ca nr) ° Lea iF Suet oy Dee oer ey i 0 J dC Utara yey a+ tae Be ett LOL 2 rea i oo Ca ee Mah cp git ak Plt) aa) 
+ te . rr) +. ar) LP fan) Cr nr re ps a. r 1 Caen we . ert. a OP ae iy ary bale we Bee eae i= Aa a dtm) ble SU ane) CO « 
Fs i Se eee » car} ‘, A a ar Oe er) are POR OR a tree p oor CONT ey Tir Pears Ge TAME SITs ea ty or aR It a ie +o OW we Bs gle ert tat Pe 
D ' Cd bap a or, fo U ‘ ' 2 eS as n eh ', ct | ’ oo. Ca) out Feb gity oi eee a Sapa oP Keltsi eg —— ° Ha NF Sat .ectn om, Cate 
. oe 1 8 iD acd ae | Lr) fouy Lie) Wan at ist ary yaa) a ee ee” PF een a a a eT at HN ee PE enue i dat i eT ie ra a taaltal 
‘ ‘ ‘8 oa) iD Ca Lu ry on of 4 Pe an _ #8 ee + Ia era Caterer Cy er ey ae a. Sar Ng bat aa OU he abated ott * 
’ t ‘ Cr) ‘ ar) Cd ar] 1a Cary ry te ote epee Fi ORs arian) ayy os for hoe, ay oe igs a der eat] idee eae Lod Neel Alt or ol 2 ae 
o v . o aL DC) y ‘ Ua ul bi Per ris ed eye luther Smet Tee at na Care Ae Lt arr ee US or pee ae ae Lay alia met it oy 
‘ oo. . ’ ‘ ' ud oon Ee tary “ape a Let dee ar era. ae Pa ae a) i sh CE SSO Daal iy i oe 3 WE Oho, old ey fen dhe 
my u 7 ‘ + CS iter Per ‘ - UE Oe ete vier r) acl LY SCs: Ale er Che Ur | rs a) Cera ir ae tal hf eee o Om 
ar] 4, ' ' ’ ry rer . F tae Ua yr et einen) Sry 4 £ ara) On nny rly Cr Sa ae a aay [ es Lr oa Tr Sr Ppt Laka ates 1d . 
. no. ° ‘ . ‘ ry LP. re | ' Cr ard Cr Leesa Ar De Meer Paar Y i) oa ry - i A aL di 3 Coa ed ar ene hs Sepp MOR NO ose eae ba 
er) e * ry ‘ i] [an rae nt | Pa ea sear +6 Oa 4 rl "@yuee it ear) ¢ t aa ry Cari rary ‘ 44, | a hae # Fage aaa Le. rs tee bend FO Om ee oy FOP mas 
i ae a 7 ' a er ‘ Cr) % © og fl choot, 6 day the on Legh tI US Ta merry or ooo ay stil I Lo) Sal etd beet at eel 
oe Uy ‘oe + Cn in or) CAA reer ‘ Oey SP ria lee ALA CBr area ie hy scr rare ee rae CP er ar eee i, i id a LE Pe eh Te te Pe 
an LU oy U n ‘ a 1 ‘ «4 a ae) "4 Char a ay CT ar “eee eee) wre, at Ls Wa ool 3 Pie, yw pe ole ek at ed htedehadal ate BL 
rs r ‘ ry ' Cee Py ry a ae) Wester ere} ry ad Cmerd Lac i eed Ch ian OR Le Pee ae Ww) +e? gas el at A Te) 4 t/a en iil lk Tae panto 
‘ mt . Ca . s rT oe ‘ UAE Baths rrr Ota a har 7S ney i ne f Ae 6 a ee T) ONO 2 OLS Fh a ga ba Mal eek dk ana 
. ) 4 D o ry Py rar) Caer Le ee or oie ier | e , 7 CC ain rien at aro vrs yt 4 ry pe gC AY et i ie ers ad ha tee ot rts rane 
o 2 ra ca of ie ny Ca ; A a o Oars cde let Tea een ON Le ty DOOR ete) Dr otic) ae Lal 2 SO ot Laake lal ea eT Penn oe Coded 
fs. ae ny i o “ : ee el Ver ear NW ag CC aD e of pee Un MOTE Ar iad Py sot WN Oe whe 5 gs ha re a O Nee oh gt eo gn Ie et we 
‘ ’ ' gies ry ces TER Tarren | ’ Py + YY # amy rar a UO se earatginensy ad SEL TFs Pee OL a ee ae ereay We 7) ON) 6 ames mye hele teed (ot ei fi 4 Ls 
1 fo 6 8 ae | u U CS a emer ta) rey oon a) a LRU) RCTS iy a, Pe combo Pe ee Ae. Wine rahi LR LET I ee Preis wn 
ea ror A 0 oo. ry m Cl er eC ‘, nr iar) eer er a a) Le he aan ’ Oe rey ery ee re ie hot ry ts mg Pl Paes ee) hel at] a ae Le 
U U LU bal ’ U oe Lar) + 6 ] * a ee oe ens er | Cae er ee er er an) Wary Peel, s om 4 en’ see EYE Tre We) ra or ere] el dahick De SIA Ve pepe 
a f . ae 0 fn . 4 A Cee rier SF Ay Ser a iar pietTEn Oe OAT ae ee ee a eC Pers te ee Or Ca eo Poth Pada a OES ore 
iu . * : o ' a Aeon oan] Te arpa) ert MOY 88s 6 te ODN og [a er er Pe Tray Cd wee ArWWotmweune wm Wh nthe eed wey 
LI 4 an + oe | a a | > para ee ICR are ae ae era ° Py 9 0 pres fox Lae Sil LC i a bade Ld Ye asy! Phd @ oem Paella obtad fk Sate 4 
‘ ' ‘ ‘ ' Ci ary Cr ry ff ny a Ae eer er So Jvae “patsy ay 44 Oi Ly Pan Pic ar dae LTT at Oe ory} Oniane ws pdfs pid eT, PL as i ert ain: 
mui ew " Leal er Pen acon yey rare nr) | ol La ot We Yh) @svaws hill Dae Ae OY had NY Ae CTT ft  petaiaedel cat LTT Perens fot 
. ‘ ) ne wak es. RT hs iu Ti ier ” i ° ry Or ee) Vou ay) haat SL re eer Se 2 a ol Le) Js eta pod me Lee Ie ey A Plead dideel af Hh 
p o " o dl aa q n n teen MEU CA WE is aay ‘at ae hed i uae ar Ve Pe RO Mio on Nad ee ee Sak ed de ete hdvghirh Lh” 
° - fn ‘ + ‘ ars ar] ‘ ra ry Ca et ee or rena ee ry Rene AV es my 2 ear ak ee Wr ee bt oe Te a ee ite wkdad oY) oh et 
i u Uy 2 a oz ee aay Pe | ey Ce eae Yo ape y \ die De Yar} i er ee 1 ee eT ee AD le od et tite Mad at hal det oa 
. ar a) i] ‘ cary r ee ee . Fy Caer ar LA Marace ae aa CROC y eae A FN Pe spew inns g Lo ae: Cr a Lode ok wey) ey Cw 
1 rare ce at ay Tae ‘ ar] road Caer] eer Co Sie arene. Late at Bart re One oe  & 04s § wr witopem » Oe po Veink gry le a) yh a 
oO iD teh eae ae NW ‘ La ON eg * if PM Seyi ia na eet OST pews wes BS acl od a ee a ae a et 
‘ 1 ‘ ae or ope 1 O ’ . Ls ere Per ir] 1 + . a nd 1 @ ar a | V a4 Ce el Pa a} Veouw 9 Lay lu fA 
ar) er 0 Cet ey 1 ‘ ar) Fier Cy A A eee oy er wrace & Hervey LN a Le eS Se alr ie bt Me at oe ts bp b) drat, ahve 
g D cp een Cae? Ih or ‘ ne Nex Le yo UP Ka eR re rrr fue s oy Ah Te SC a ot MT ‘se i lla he ee tte hbk deol! We YT) (Mise ciated ol Ah al 
an) . ve CS eery lee] LS) ie a r o a "done ueaea, Yow Ay a ee Pe TT Ct Perera et a a Le YE heed OR TTL Say bot wade f beth Died fo ie ol ol) nd 
‘ ar ‘ bet a Mee oat fe 7 a 0 ard ‘ A fom yp ten ee AC Day eo ae Ont ene a ae ee meee y ee er Areal or etn | de an sate be 
. ian) . Md U Cae | ‘ Ca rr a ? Lie rcs CLC rere ae Soop gy itt a oa PN OR oe oe pes ‘Porias (ett PAAR Sieh eas d sd 
D U ' ee ee RS ee ee ee a ae Oat aCe ee * Ce Ne we ey ty ego TG eee wre faeces ge hes bn Mi padre 
A U Ma U ae i; wa 1 Py ‘ A hoe er ar er ery ar) te ee ee i ee et “Pee ryegs, etek a ok} i idee ST ce een eee 
7 4 U ea i ae Lr ar) + (es 66d hte vue ' CT LC ler ec # Rr yey Che Le Tes VHP ws ld bie hue Meh Naat Mh eh alee oer ys 
‘ Ua | A i] fa UO Le er ee Yee + ae oy Par Ca Tee emer eee ra ai 6 Oe A tl Se Pe | Le ity aa ke i epny ted = $ Fem im pen 
p fl p ere q f Ae a Pe orn Ciacci ar ie 4 nen Tarr er eer rer Nd PAL Peed Lat ei eaten, ere wt 
o 4 ane f Ca . ar a) Te aeey Pia eed lar Sasa a f° Ok eee oso ‘ Tae nea BIO RO rs Fata 
me ' ar o ae oO ‘ r _ err) ar ey Cn a Te an) mM Ti Or be sete VV ww bee a LL, are) Rate dee Eh ee, BT b 
so o y ‘ fas Cie JL a bd pmo scat et haat it Jar ian Le eer Sepia i Y betty» ee es OND ony Ae glen Ne fem oe 1 we + We gt w 
o eee 2 H tana ay : 0 r rr LU Coes 04 Behe gen ay ee ee ake SIPs hE ih aed aw YY on 1) 
Fy ’ ' Cn a] . rar + oo, a] t rin] ey Meri ery ‘ ry vy a Ee ee oe eS cae oe ny) OO REA RT Lene eT ida lag’ Mk a) 
‘ ’ u eae naka) eee Le cD Pee a foc) 4 Ca | P of ta 6D bene ‘a+ © Pt Aa tt ae aie ROT) aN or 
Cie ane ao o iD La iris bear Wit ale ta oY erate Oka MOR ae ely Cee YA dd eer ee ee dl aR Sk ed 5 
1 A 7 . ar A ae fry ' Ue ree ar | er Lt a ry ee a er i ee ey (Fi Few sn ge i eb ee a) bd aera ACS Tee St all sulted dcheol 
‘ ‘ " a? Oo | ' iu ' ri oO a TY Oo * OIE LIE BPN OL ir) ofa Lh et ee Se CO Sate ey ceded 
D 2 A ‘ cd A r e r ae ay A ReCeCU a, eeamriey A soar Corie uty We agg Cat a en Pt fete ee“ Pvpwe, > adie) WH Pehl) re] Nett Dh SY a oh al Pa 
4 al a t) ’ a a ey ar] Pa i} ta ’ Ces Par ihe : Ce er . 2 a te) © Cot wh ges, o@we Lal Ae! Laan rf aA dd Cd Ae od ae oe ee ol 
Car ’ ta . oar] ar . ar] err, ry ‘ A aT Sar om) Cn a ry age va Par cera Oe memes fre 88 EA IATe fois ap Nd ahd, ad te are! 
rn r ’ en A o. f Car Ce a ay ’ rs F A + ree Ce PL Sd an en it a ey We ne Lier eer wotwe etialaak ot ee vive « 7 oe ke 2) belt A a 
o ‘ eee ee ie ‘ eo at Ar ary f Ue ni eRe he ee Te ey FA ae so me i CECE elie Pee SA boldhil a deal 
4 +e 8 o8 see a ” Fi a nC a Pe a ‘es ea wy ar Le ee a Ae Oe FY 0:0 9 regs © aN Heat geacy tga Laas Lilet tly nh 
A o so. sf soe Lares -2 he r eS Hien % toy Ch Ee eA Se PS fe bh OND» wHE,y IPE Pe ed Label hots ke 
i et mON au en BLU We a ora "e494 re es ary ie ee a une haat PI kd BOON LAO I Ut as 5) EC WA Oe Fe een S Sethe ad Aad, 
a E ‘ : U hue tess. ol es UIDs Train IAs ee CE aca r) a ee 86 tes? Pew we Utd PrP sy © rw bie et ft at hat liad he 
vo ' Ud el ee CLs eer > be a i On eee ettarith el Dire eee a ol ey a) I ad ade A Peal os AA Sl a on SN riphielegtiah Mak | at 
rr r ss hs eet are | s t ele le Oe Yang | 4 ~ of 8 a | eC A a de i Doar ee pet rag Gy age: Sed MMe add ha 
rt pl Aes. areca 6 0 6 PRU ie bas fa hah, via eg : ee Kr tae of ry rn ae | Ch ee PPT OWN my PERNA 6 wwe at Sane A, ae, 
ager . O ' er i Let Me oe ta Uae tina r n) Cr ne rhs 0 Tore er pa Fee ay a Lee 7k aes bd PA eS te Lf Od, a ied dood od 
o ee) fl fl aL areata eA eT AS a oan Cr) a  . P ) Le i es " Saugpray eA ak a ie | CAL ht hee rs wae | ahd dod Cp deplete he ahs. Pe 
V+ z : nee ‘ ate oye oer mer Ayr te es Led | TY saan a A at PO ater eR CSR Pee. hol a erie toad 
MN ne Poe U 1 ar ae ‘ om ti rn’) ar) ower fn) Cae ey a) f Pee pira, tee TaN Sedge | A Ce | Cd ie 2 ee iweb ee) Leh od LT Flop he biel olla ok 
oo ‘ D fo ber Wer Ter ot fy | ’ fog ie ea) ro ao ov. me r Le | feng Cee 2 we) Ads ee) ea tell Avot he fo ‘YF oc bepy 
D a oa oar fa ts oe 5 eee le Ae heey. ny ee Si Neer aed MD oe 2 BON EERO Fare oe yy, 5 
D ‘ ‘ Feet enon sos ay a ee er Le ahr 2 Te we) “hh beRe S48 Here ated Bs ary) ieee STS ey ere 
sae a D t Pt a H CI Fg ON 6 + py ; ee ee ee ood +5 eb im Yemen a, POV OP a EO 5.5 ee ok Ms tr) i 
eae 0 0 $ i n ‘ Le ra o Oe Co et 2 a ay Be Ih ae ee ee Ea ae Teer ee SOE Pe Allele dag) a 
soars , noone a) meer) D Cee hee F t. cena O Le aed A, 2 rer eg i ah ea LS oy Yn ot Mh So Tk Pra a att aeharstidl al ot af 
A ‘ of Un eae ‘ i ‘an ‘ ae | r io te , ¥ Le eer Pere a CSE ie ¥ Oe de ae ic we ts bh Ded ol ee ee | Sl a te es a Ot 
ee eee is Te er ria 0 ree SUS f ce ee ier ) a ee Jer ae ei oe TTS WT hd ae PT | erg t bel, nflinclia Regt A Le os 
A o D iy P a hata US exc oy ® , re | iia rei OF We ie wy ent fey wy ‘eas te ee Tht ™ iL. 
0 ei ' an ieee ee i ee Soe tliar a t ae A ae ee ee Lit hE Tey pad Oi RRL Ripe yg 
‘ 7 y to 4 Cty arian arity yeaa ‘op Poty 04% weer ive, Meee del ee Hi hebt F Md ut Le cad pt rapes Leah poe 
f iar a ae oes F A rr way és Date CT ar sy a a a Zs bh Ao oh Ae of PRETO ES py Lie oO AVS WF By rng: 
, , "4 an o oy oo Ste n ‘4 To 6 Oe pe My vie >, er er Fe LP Le hat os OE erty A dtdmth ad at of, 
A D 5 On ra er Py) A ae te the A n 14 Ve Ree ET ¥ yee tg hn) hse PRS PP srt wy 0 Lhe weve 
MN ‘ A Ps a ae Fi ae ; a Ps et ee | aes + eh Mets a Ae oe} bd at a oe S@e-. , Pee 
g eee oe Pree Paar my dey hype (ope SOR AS Faso g ed we ee OF ea 
. . o ’ eee F ’ F oy ee Poe r a 0 ar ew a he er ee eer ee 
) r ) Fy r ‘ ay Dear Stes 
- Car “oun (Oa Decuria OPIS RT ver 
* e 
ery ter 
4 oe 


DU a> | ee 
xs — -_ 
e ' 4 omy & P84 

















NAVAL P 


Monterey, Galiiornia 





OSTGRADUATE SGnUDL 





Ht baht 9 


NORMAL APPROXIMATION FOR 
RESPONSE TIME IN A 
PROCESSOR~SHARED COMPUTER SYSTEM MODEL 


by 
Suriya Pornsurlya 


March 1984 





Thesis advisor: Be Jacobs 


Approved for public release; distribution unlimited 


1215671 


DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 93943 





as 


Sas = ae - 





SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


ae “et 
READ INSTRUCTIONS | 


4. TITLE (and Subtitie) 5S. TYPE OF REPORT & PERIOD COVERED 


Normal Approximation for Response Time in a Master's Thesis 
Processor-Shared Computer System Model March 1984 


&. PERFORMING ORG. REPORT NUMBER 


















» AUTHOR(e) 8. CONTRACT OR GRANT NUMBER(2) 





Suriya Pornsuriya 










| 10. PROGRAM ELEMENT, PROJECT, TASK 


. PERFORMING ORGANIZATION NAME ANO ADORESS 
AREA & WORK UNIT NUMBERS 


Naval Postgraduate School 
Monterey, CA 93943 










12. REPORT DATE 


March 1984 


13. ae OF PAGES 





- CONTROLLING OFFICE NAME ANO AOORESS 


Naval Postgraduate School 
Monterey, CA 93943 







1$. SECURITY CLASS. (of thia report) 





- MONITORING AGENCY NAME & ADDRESS (if diiferent from Controiling Olfice) 









DECLASS 
SCHEDUL 


1Sa. 







IFICATION/ DOWNGRADING 
E 





- DISTRIBUTION STATEMENT (of thie Report) 


Approved for Public release; distribution unlimited 






17. OISTRIBUTION STATEMENT (of the abetract entered in Block 20, if different from Report) 


18. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue on reverae aide if neceeasary and identify by biock number) 


Response Time, Central Limit Theorem, Heavy Traffic 


20. ABSTRACT (Continue on reverae eide if neceeeary and identity by biock number) 


In a time-shared computer system, the processor allocates its processing 
time equally to all jobs submitted for service from a fixed number of 
terminals. Under Markov assumptions, i.e., independent identically distri- 


buted exponential terminal think times and job requested service times, the 
distributionof response time of a tagged job theoretically can be determined 
by solving a system of differential equations derived for each initial system 
State. However, explicit closed form solutions to these equations are quite 





DD van 73 1473 EOITION OF 1 Nov6S 1S OBSOLETE Unclassified 


S/N 0102+ LF 014-6601 1 secdRity CLASSIFICATION OF THIS PAGE (When Date Bntorec 





Unclassified 


a SS 
SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


complex. The Central Limit Theorem and heavy traffic arguments suggest 
normal approximations to the distribution of the response time. Simulation 
of the response time is used to study the accuracy of these normal approxi- 
mations to the response time distribution via moments and quantiles. 
Finally, the analysis is extended to a model for a system with two types 
oieeerminals. 


S’N 0102- LF- 014-6601 


2 


SS EEE Sa 
SECURITY CLASSIFICATION OF THIS PAGE(When Data Entered) 








PWeeremeG@meor DUbDiic release; distribution unlimited. 


Normal Approximation for 
Response Time 
in a Processcr-Shared Computer System Model 


by 


suriya Peornmsuriya 
Ensign, Royal Thai NAvy 


Submitted in pa 


Tate 
requireménts 


ye 
C 


a u 
mor tc 


MASTER OF SCIENCE IN OPERATIONS RESEARCH 
fron «he 


NAVAL POSTGRADUATE SCHOOL 
March 1984 





ABSTSBACT 


In a time-sharec computer system, the processcr allo- 
cates its processing time equally to all jobs submitted for 
Ez Markov 


b 
service frem a fixed number of tsrminals. Und 
assumpticns, ise. independent identically distrib 
nentiai terminal ‘think times and job requeste 
times, the distribution of response time of a tagge 
theoretically can be determined by solving a st 
differential equaticns derived for each initi 
state. However, explicit closed form solutions to 
equations are quite ccmplex. The Central Limit Theorem and 
heavy traffic arguments suggest normal approximations tc the 
distribution of the respcense time. Simulation of the 
respense time is used to study the accuracy of these nertmal 
approximaticns to the response time distribution via moments 
and quantiles. Finally, she analysis is extended to a model 


for a system with twe types of terminals. 


| * 





DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 935943 


pape Of CONTENTS 


1 INTRODUCT ION e e e we @ e @ @ @ @ @ @ e e @ @ @ e a 
Ae EACKGROUNC @ a e e e @ e 2 e e e e e e @ e e @ 7 
Peete NG SYSTEMS .« « « «© © © © es wo wo 7? 
Core eb ee RNDeee ee ROX IMATM DISTRIBUTION . «© « « «© « YQ 

Ii. ore ermeieoerOn A SINGLE JOB TYPE MODEL . .. 12 


RP ROIUTGTCNN) fe ce ew we ew wt we ew 42 
B. CONDITION CN REQUIRED TIME AND SYSTEM | 

PE es ee wt wt tw wl es 13 
C. CONDITION CN REQUIRED TIME. .......26. 14 
D. MOMENTS AND VARIANCE OF RESPONSE TIME .... 16 
ROOUNUMESRIRETNURESULTS . . 2 . «+. « « « « « « « « 18 


i. See toro lmOrmmO Nn OG TYPE BMODE Lis. « « «© © © oe ew 20 
A. eet UGrlmNe cums Seilisialck 6 sulisuleslls « « « e « “20 

een CeCe er CCRC Uc cae 6 = «© «@ « « © ©« © « « « 21 
Cont DONE ORAS 2Z-TEBRMEINAL SYSTEM . . « © « « 22 
DPERRAECORINGUIMGE Ts 0% 6 « © © © © © # «© s © © « « 22 
Pemeeibemllaa: RecultS . 2. « « « © « « © « « « 24 

De SeeuwetnON PORTAN N@=DERMINAL SYSTEM . « « « e 20 
oo G EEC gGGWNsIs =. 0NRN «06 0 «© « e 29 


ZO M ete mO tenes PCNSewraMe . « + « « ©« « « « 28 
Se COmputation Cf Simulated Variance .... 30 
4. Computation cf Simulated Skewnesess ... . 30 
DP ecommeateontnor simulated Kurtosis .... 31 
6. Standard Errer of the Mean Response 


Time @ & @ @ a @ @ e @ é & e @ 2 Cd é @ @ @ 32 
Vemeppitlem=ical ReEesSsultS 1. « © « © « « « e « e e 33 





aaNe 


VI. 


tener nO a ATION FOR RESPONSE TIME . « « -« 
A. Wee eh ON = 6 © « «© «© « « « # © es « « e 
B. APPROXIMATION BY CENTRAL LIMIT THEOREM .. 
i veeectieso lames Theorem for iW¢t) . . -« 

2. A Central Limit Theorem for response 
Pewee) so bes s « «© « «6 6 © 6 6 le 
Coe Oo eee TGN BY HEAVY TRAFFIC ANALYSIS . 
De COteeereesc Ooh TO SLUULATION DATA . «ce « « e 


MODED@ren er si Stites PIn Two TYPES OF TERSINALS 
A. re UNGrPMUOTIN Ss) 6 Gs 0 «© © 2©« «© © «© «© @ «@ @ 
Eee cote DISTREIEUTION «© « © « « © @& « 
Co Amer LEM THEOREM FOR THE RESPONSE 
PL OSPctisiNs Nellis 2 e ee « « « « « e« e e « « e 
ID) EEAVY TRAFFIC APFROXIMNATION FOR THE 
PemSsP Ger thE .« « «© © © © © © « © © © © © © 2 
i. Sew HeeenON GMs «0 © « «© © © #® «© © © «© «© © « 
(PeeCOPEMLe nN . SMS 5 ee « eg Ge Se lw le 
Ze Computation cf Simulated Mean, 
Variance, Skewness and Kurtosis ... 
Seeeewiewemcal ROSMItTS . . 6 « «© © «@ © «© « 


SUMM ARY AND CCNCLUSI ON Q @ @ @ @ e e e e e s @ 
A. SUMMARY @ @ 3s @ © @ ® ® @ & @ Q e o ® ® a 
Bus CONCLUSION S @ e 2 e ® @ @ @ ® e @ @ e @ a 


APPENEIX As SIMULATICN PROGRAM FOR ONE-TYPE MODEL . 


APPENDI 


Bist, OF 


INITIAL 


X Bs SIMULATICN PROGRAM FOR TWO-TYPE MODEL. 


REFERENCES 2s ® ® @ D ® @ og 2 s @ ® ® ® e ® 2 


poorer WON Tiers ss s ss 6 elt lel tll lt 


56 
SO 
5 
oy 


39 
40) 
42 


50 
50 
51 


wD 


58 
63 
oo 


65 


68 


9 
1s 
8 1 


83 


se 


104 


106 





I. INTRODUCTION 


A. EACKCROUND 


Pes 2ertinel, “submits e job 


3 


mien a computer user, f£re 
to the ccmputer, it would Le desirable for him that the job 
ke precessed right away. However, this rarely happens in the 
real world, because a computer system usually consists of 
only cne Central Processing Unit (CPU) and many terminals. 
So if more than one féerson uses the computer, there will be 
jobs that request precessing at the sama time. Hence, scme 
kind cf gueusing system or time-sharing technique will be 
needed +o crganize the allocation of processing time to 
those submitted jobs. 

Generally, a computer system has twe types of precessing 
time allccation. One is called "Batch Processing". Ali the 
jobs sukmitted this way forma kind of queue and wait to be 
served (processed) according eo the well knecwn 
First-Ccme-First-Served (FCFS) POLLEY, 1. ¢. the first 
Submitted job is precessed to completion, then the computer 
is set up for the next waiting job; in some systems a 
Eriority policy based on the length of a job is used to 
determine the next job to get dedicated use of the 
processor. The other one which is of interest for this 
thesis is celled "Time-Sharing Processing". Tt is based on 
a technicue that permits concurrent processing of twe or 
more jobs. Each job no matter when it is submitted gets an 
equal share of processing time until completion. 


EF. FROCESSOR-SHARING SYSTES. 


The so-called "Processor-Sharing" or “pure time-sharing" 


in computer engineering is the system in which the processor 





Shares its service (rrocessing time) egually among all jcbs 
submitted. In other words, if an individual job reguiring a 
certain amount of precessing time 1s tagged and submited to 


the system and finds (j-1) other jobs being ptocessed 


= 
ct 
is 
(D 
is 


tam new cn all j jobs, Will each receive service 


(processing time) egual to (1/j)-th of a (processing) <¢ 


c+ 
c+ 


ed 


jobs receive service changes each time a new arrival jeins 


Dame fer time unit. Gfecewese the rat2 at which submi 


the system and each tire a completed job departs. 


mies acceadeemoneG: CONCUTEL Capacity allccation may be 
described in more forral terms as follows : 1£ the chance 
that any single job, processed alon2, finishes in time 
interval (t,t+th) is eb + oh), ( €xponential-Markov 


service ), then the chance that @ particular "taqged" jeb in 
the ccmpany of (3-1) others finishes in (t,tth) is 
M(h/ 3) +0 (4) as h approaches Zero. 

Clearly, there is no waiting line in precessor sharing 
c£ the abcve type. It permits short jobs access to 
processing right away even if they arrive after longer jobs. 

Precessor sharing is an approximation to the processor 
sharine “Round-Robin" model. In this model once a particular 
tagged jcb enters the system, it joins the end of an ordsred 
gueuve. khen it reaches the service point it is allocated a 
fixed quantum (q) cf service time. If the job completes 
within this time it simply leaves the system. If afterg 
seconds it still requires more service, it is immediately 
returned to the end of the queue. This process then gces on 
and cn until the required service is completed. 

In tke limit, however, aS q approaches zero, the 
Round-Robin system becomes the Processor-Sharing system. The 
latter system here will be studied for its characteristics 
reélatirg to the distribution cf response time, ise. the time 
that a tagged jcb requiring a certain amount of processing 
time actually taken tc complete the service and leave the 
Systen. 





C. MCCDEL AND APPROXIBATE DISTRIBUTION 


Peete wine 8tttrs. Study of dslays to arriving jcks 
Deets =Oces=o" Sharing waS conducted by Coffman, Muntz, and 
feet crer: (1970). [ Ref. 1] They assumed a steady state M/M/T 
system under processcer sharing, foo. POSS eOnwatravels and 
exponential sarvice times with a singl2 server processor. 
They then derived an expressicn for the Laplace transform of 
the waiting time distribution of an arriving job conditioned 
en the processing time it requires and the number cf jobs it 
finds in the system cn arrival. 

The preperties cf the response time, Ry! waeven’ )2he 
processing time required by the arriving job under processor 
sharing system was further analyzed by D. Mitra (1981) 
(Ref. 2] for the following model. A system consisting of N 
terminals and a single computer (CPU) can be modeled as a 
classical machine-rerair situation : each thinking terminal 
(failure-prene machine) applies for computer service at rate 
\ - and queued or waiting jobs are served at rate pH as 
long as any jobs are present. If Markov assumptions are made 
throughout, then X(t), the number of jobs at the service 


stage, ig a birth and death process with transition rates : 


R(t) = J---> X(tehy = Gt Tos Ah + oC) 
“==> K(teh) = F- 1 oF pyh # ofh) (1.1) 
---> X(tth) = 3 ee fre + 0 (h) 
where A; = A(N-3) H= the rate et which a job is 


submitted tc the comruter when there are already j jokts in 


the system, and Bi = H for j 2 1, otherwise being zero, is 





the rate at which tke computer gives service to all j jebs 
submitted. Based on the above transition rates, the distri- 
baeLon cf response time under processor Sharing is 
characterized, and the mcments such as nean and varience are 
Emma Under iwnterestang conditions, i.é. themCcOndi t2Ona ! 
respcense time, given only the precessing requirement T *ime 
Ua tS!, derived from the condition that the tagged job 
Earring T time Uunats of processing arrives to find (j-1) 
cthers in the system. [Ref. 2} | 

Gaver, Jacobs and Latouch: [{Ref. 3] have generalizec and 
extended the previous analysis by introducing the idéa of 
processor Sharing in an arbitrary birth and death precess 
envircnment, thus allowing quite general t¢erminal-computer 
interactions to te represented. In the process, the meaning 
cf "system state at the moment of tagged job arrival" is 
also clarified by Lavenberg and Reiser. f{Ref. 4] Response 
time characteristics are computed under the assumptions that 
rrocesscr-sharing service rates are processcr-state- 
dependent ina mcre general way than that described earlier; 
this allows for approximate representation of overhead 
Fenalti2zs and also of job scheduling. Other characteristics 
of tageced jcb response are also studied, e.g. the accunu- 
lated precessing work, W(t'), actually performed on that Job 
by elapsed time t*, with t' < T (T = required processing 
time) fcllowing job introduction; note that W(R) = T, So 
the first passage of W(t') +0 T ais actually the response 
time. Although differential equations may be obtained for 
transforms cf W(t) under various initial conditions, and 
hence, implicitly for its distribution, the results are far 
from being ¢xplicit and informative. However, central limit 
theorems for additive functionals of Markov processes, or 
for cumulative processes, allcw the conclusion that the 
accumulated work acccmplished by fixed time <t' ona "long" 


job is also approximately normally distributed. 
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Mae etencdiy, < ZOLMal approximation is shown to &@ vais 
for the simple model --and probably for others as well-- 
uyen the number cf ccmpecing terminalis becomes large , i.#%. 
Pie teneavy trattic¢ conditions. The quality of 
Bupeeminao ens fOr finite job lengths and for 
number of tsrminels will be assessed by cimulation methods 
in chapter IV of this thesis. 

The differential equations for the méan and mements of 
the response time of a tagged job reguiring a fixed amount 
See neecrpesing tine, given thaz it enters to find an initial 
number cf cther jobs being frocessed, will be derived in 
chapter II. To remove the condition of initial system state, 
we Will use the steady-state distribution of the number of 
jobs at the service stage. This will also be explained in 
the same chapter. A frocedure to simulate the response time, 
given again an initial system state, will be described in 
Chapter IIl. The empirical response times obtained from 
Simulaticn will then be considered as stratified randon 
Samples. Fence, a method of computing the central moments 
will ke given accordingly. | 

In chapter V, we will study a bivariate birth and death 
frocess model for acomputer system having two types of 
terminals. This model allows relaxation of the independent 
identically distributed exponential service requirement and 
terminal think times of the model described by Mitra. 
{Ref. 2] Under the same conditions as the previous simple 
model with cne job type, we will derive normal approxina- 
tions for the distrikution cf the response time of a tagged 
job reguiring + units cf processing time. As before, sinmu- 
laticn metheds will b2 used to assess the accuracy of normal 


approximaticn to the distribution of the response time. 
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‘oe 7 
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i i ' 


II. NEA 


AND SQHENTS FOR A SINGLE JOB TYPE MODEL 


A. INTRCDUCTION 


mets =) Ciaeeet, we consider the birth-death process 
model with rates in (alee |e Ever though ta= unconditional 
mean and variance of response time of a Processor-Sharing 
system nay ke obtained by deriving the Laplace transform of 
Premedudiictrium waiting time distribution fRef. 1] it is 
also interesting to develop scme results for the conditional 
expectaticn of response time of a job requiring T units of 
processing time, since in the real world one might rather 
wonder hew long a job that requires certain amounts of 
Beoecessing tmme will Ere delayed after being submitted tc the 
systen. 

We will show that the equilibrium mean waiting time in 
the FErocessor-Sharing system varies linearly with the 
service time requirement T, i.e. E(X) = AT/{ lA (1- A/ yr: 
as T ---> ©. Thus for arrivals having a service ‘time 
requirement less than the average, i.e. T <K 1/M , the mean 
response time is less in the Processor-Sharing system than 
in the First-Come-First-Served systen. 

Ic derive the ccnditioral mean response time, a given 
tagged jeb, i.e. a particular job that enters to find (j-1) 
ctherse waiting for service time, and that requires "T" units 
cf precessing time will be considered. Then under Markov 
assumreticnse a system of differential equations will bes 
established to allow computation of the conditional mean 
response time by numerical methods. Other moments nay also 
ke cktained by mé€ans of solving other systems of 


differential equations. 
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Eos CompitztOnN OW REQUIRED TISE-AND SYSTEM STATE. 


a 


Under the assumpticns of a Markov process, i.¢. [rirt 


and death process, on vthe number of jobs at the servic 


rt) 


Pregence @ Ssys-em ce "NY terminals anda single ccmputer 
(CPU) with transiticn rates es in equation (1.1), a system 
cf differential equations for the mean response time may be 


derived as follcws. 


Let R refer te the resrfons2 time of a newly arrived 
job, and 
pe = ECRIX(O) = j, WOR) = TJ], (2.1) 


fie cenatticnal empectation® of the response time, given that 
the tagged job is initially in the company of (j-1) cthercs, 
ete arrives to find (j- 1} jobs present, and requires 
“work cr processing time equal to fT. 

Let dj and Mj be as in (1.1). Consider all the 
fossitle system chances in (0,h), and subseguently; any of 
the fcllcwing mutually exclusively events may occur: 

(a) new job arrival, bringing the state to j+l, an 
event cf prebability A; he 

(tk) accompanying job departure and return to think node, 
an event of probability (J- 1) BR (c(F)7HD 4b. 

(c) no change in accompanying system state but a reduc- 
tion in remaining taaged-jok service of (r(j)/j)h , an event 
me mecbart lity 1-( rj + (5-1) ul Ca aye) hs 

Ali cther possible events are of probability o(h) = and 
may Ee icnored. 

The term r(j) used above represents the fraction of time 
the processor actually spends processing when there are j 
jobs being processed. The fact that r(j) < 1 £4xrepresents 
overhead. 


tok 
(a) 





wetting U, 3 ea 7 3) (Meet opmand (cy Lead tc 


m,(T) = hem, (T-(r(9)/H)B VEIN“ (AY * By) BI (2.2) 
+ Ayhos, Crs (CS) / Hb) + Myhay (T= (279) h)+o(h). 


Subtract m4 (T- (r (3) 73) 4) Peomecach Side , ther divide by h 


and let h ---> 0 to get the differential equations : 


(HS) 7Hm(T) = 1 - CA * PAD BD (2.3) 


; 4 
+ r,; m4,(T) + Ap my_\(T) - 


This is a standard system of linear differential éequa- 
meon=s , inatial conditions are m, (0) =O fOr ala -4 x. A 
solution can be cbtained in terms of Laplace transforms, by 
exponential formulas involving matrices, or, numerically, by 
use of standard computer ccedes for the solution of systems 
ef linear differential equations. Sanples  explzcit and 


comprehensiktle closed form results do not seem attainable. 


C. CCNDITICN ON REQUIRED TIME 


The condition that the tagged job entered to find (j-1) 
cthers in the systen, Bie anes X (0) = j, can be removed 
according to the stationary distribution that corresponds to 
the system state found by the arriving job. The resulting 
expressicn allows the conclusion that the expected responses 
time is "linear" in the required processing time, T. The 
result here holds for quite general birth-and-death precess 
model, not cnly for the simple machine-repair setup detailed 
here. [Ref. 5] 
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The derivation cf linearity of the expected tresronss 
time is developed, in outline, as follows. 

First, observe that the long-run distribution of X(0), 
@. the number of jobs present (including the tagaed jcb) 


just after the teggec job enters , is 


ee CTT; Ajey : e Wer (9) ¢ J = 1e2ee0-N (2.4) 


where c is selected sec that the q 's sum to one, and for all 
ws, wT, = THe (AoA = Ajay) 7 (HM, © = BY) Wome em alone: y 
weoummbpucton (essumed to exist) of the Markov chain X(t), 
i.e. the number of joks at the service staqe, with rates as 
in (1.1) with Mi = HxI()- The equation (2.4) is intui- 
tively apparent, fopmcimee es rond-run probability that a4 
Beeemeiticn from j-1 t¢ j occurs in (t,tth) is TT; Ago as 
h ---> 0 and hence equation (2.4) tolleows by nermalization. 
A formal preoft can be provided based either upon an embedded 
Markov chain formulation, or upon the theory of additive 
functionals of a Markev process. Weert. Ojone distri bucion 
$q,$ has also keen civen ty Kelly. [(Ref. 7] 

Next, use equaticn (2.4) to remove the condition that 
aoe = 3. Put 


N 
1) =E Ef RIX(0) ,W(R) = = partie. ‘ : 
m(T) = EY ECRIX(0) oW(R) =T] 2930 (2.5) 


Then in terms of the differential equations (2.3); after 
Soe reengueneough ky j/r (3), one obtains, Ween 2naeral 


conditions m; (0) = 0, 
N N a 
m'(T) = (j/r (j se -(j/T (jj -(\:+M)ns(T 
= G/t (5)) 4, 245 3/t (5) 0 - (Ag+ By) 25 (7) 


nS | 
+ ema (D * My my (T) I 


KS 





N 
Rea) = es 725) 9; - (2.6) 
J2A 


mame 3: fellows that the long-run conditional sxpeccted 
Beapoense tite of the processcr-sharing system is linear in 


the processing time requirement T: 


NI 
EEE) =") => T3720 (4)) 4; 
jos 


TE( X (0) /r (X (0) ) J. (2.7) 


Apeerently no such simple form exists for Var{ RIW(R)=T], 
althcugh Mitra [ Ref. 2] has given a formula for a particular 
case. I+ will be shcewn, however, that the above variancs is 


mrdemi@mproporcasonal tc T if T is large. 


D. MCMENTS AND VARIANCE OF RESPCNSE TIME 


The ccenditional moments fer response time of a job 
requiring T units of processing time may be obtained by a 
Similar derivation te that used for the expected conditional 
respcns2 time. For example, to find an expression in differ- 
ential equation form of the second moment, one has to 
ccnsider all the possible system changes during time period 
(oh) as has keen done for th2 mean response time in 
rrevicus séction. If the conditional second moment of 
respense time of a tagged job that requires T units of 


processing time is: 


m2, (T) = ECR21X(0)=j,W (R) =T], (2.8) 


then, the fcllowing results subsequently occur: 
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m2, (T) -m2, (T- (2 (5) 73) h) (2..9) 


ad 
po eee eae) 3) eal Ay t Pg) 27; (T= (2G) 79) 2) 


iS 
+ AymAy (I- (J) 73) h) + oa Gaye) 


As h ---> 0 one cktains the differantial equations: 


ed 
Sane fame (i 7or) = aa (T)- (2X, + ean 
pd 
+ rj ere (T) + 5 m2 _4(7) ; (2105 


Again, the distribution found in (2.4) can be used to 
remove the condition that the job enters to find j-1 cthetrs 
in the system. 


The variance may then be computed by usual formula, i.e. 
Waseem ee (T) | - [EL R(T) ]¢. 


Likewise, the differential equations for the third and 


‘iD 


fourth mcments, m3, (T) and m*(T), may be obtained ky th 
procedure used to evaluate the first and second moments. The 


expressicns for these moments are as follow: 


(eS) 7) Cams (1) /AT] = 30% (7) - C+ PMT) 


pre (7) + ms _ (7) , (Dental) 


(x (G73) Came (1) 747] = 4m, (7) - CXS +H) 4 (7) 


wo 
+ Ayre (2) + yim (1) . (2eA2) 
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Cnc¢ the conditicn of initial system stat@2 is reroved 
uSing the expression cf (2.4) we can compute the third and 
fourth central moments by expanding the powers in crdéer to 
calculate them in terms of moments around the octigin 
cbtained frcem solving the system of differential equations 
descrired atove. The central mcments may be =xpressed in the 


follcwing forms: 


mek (1)-EL R(T) }) 7] 


EC R3 (T) 1-3E( R2 (T) JECR(T) ] 


+2( a(R (T) JJ? 


EC (R(T)-ECR(T) D4) = ELR* (T) J-4EC R3(T) JECR(T) J 


Ome een (.)) ) jes BE R<{T) })%- 


The skewness and kurtcsis of response time are then computed 


as fcllovyws: 


Skewness EC (R(T) EL R(T) J) 3 J/(Var[{ 2(T) J) ¥, 


Kurtosis = {EC (R(T)-ECR(T) J) *J/((Var( R(T) }) 2)§ -3. 


E, NUMERICAL RESULTS 


The conditional expected response time m;(T) and the 
conditional moments DAD) mn, (T) and n*,(T) can be computed 
by solving the differential equations (2.10)-(2.12) using 
numerical methods, ees €ither linear or Runge-Kutta 
rethcds. Once these results are obtained, witemc Cnid acon 
that the tagged job enters to find j-1 others in system is 
removed. Hence we will cbtain the mean and the saécond, 
third and fourth moments of response time of a job that 
requires T units of processing time. These values can tken 
ke used to compute tke central moments, and eventually allow 
us tc determine the variance, skewness and kurtosis for the 


distribution of respense time. 
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Tab I shows tke means, VeteanGc=s, Skewnesses enc 


lé 
kurtesises cf response time cf @ job requiring T units of 


system cf 2 t2arminals with arrival 


u = 2. We 


processing time for a 


rn = 1 


required work 


rate and service rate Sc aoe 


ertemececanes | 2a=gqe the distribution of the 


response time is modérately close to a symmetric one, since 


the skewnesses are small. The kurtosis is also approaching 


zero as the required work time becomes large. This suggests 
that thea dastrspucson of “the response time may bs 
Pop=oximated by a normal distribution for larg2 required 
work time. 

—_— = —_ —— i 
| TABLE I 
Numerical Results for 2 Terminals | 

! 
nN = 1, = 2 

| | 
| Zine Nean Variance Skewness Kurtosis | 
| | 
| 0.1000 Oo 13se 0.0020 JeoGs.s = 1.34 15 
| 0.2000 0. 2668 0.0074 0.7078 -1.3042 | 
G2o'000 0.4002 Os OmiS2 0:.2740153 =1.2428 ! 

| 0.5000 0. 6668 0.0358 0.6876 -1.0823 | 
| | 
| 1.0000 1. 3885 of, 1021 0.5729 -0.6024 | 
| 2.0000 2.6719 0572656 0.5118 -0.4626 | 
| B20 000 4.0084 0.4569 0.4614 -0.3802 | 
ee 


A 





III. SIMULATION FOR ONE JOB TYPE MODEL 


Aw INTRCDUCTION 


The numerical method of computing the moments cf the 
conditional response time of a tagged job that requires scme 
fixed amcunt of processing time indicated in the previcus 
chapter is generally sufficiently accurate, especially if 
carried cut by the Runge-Kutta method. However, the distri- 
bution of the response time is also of interest. We will use 
Semulaticon to study this distribution. We will describe a 
Simulaticn reutine fer response time of a jceb requiring a 
fexed amcunt of processing time for the model with one job 
type descrited in charters I and ITI. 

The ccnditional response time of atagged jcb that 
enters tc find j-1 cthers initially present in the systen, 
and requires T units of processing tim2 can be simulated for 
the mcedel if the jek submission rate to the processor of 
each terminal and the processing rate of the processor, i.e. 
~ and MW. are known. Under Markov assumptions,the number of 
jobs, including the tagged one, in the Processor-Sharing 
system which consists of one processor and N ‘terminals is 
considered as a birth and death process with transition 
trates A, = A(N-j) and B= PCG-MxrCH) 75, where r(j) is 
defined as ~¢ fraction of time the processor actually spends 
processing when there are j jobs being processed including 
the tagged job. Thus the interarrival time and the departure 
time (work completed) of the jobs in the system are exponen- 
tially distributed with parameters A; and wy, respectively. 
We use the LLRANDCMIIT package available for the Naval 
Postgraduate School computer system +o generate the two 
exponential times with rates dj and BS respectively. 
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The 


ene 
and the sojcur 


= 


WZ) 


,erated arrival and departure times are ccemrarec 
n 


time in state is determined as well as the 
a 


\-4+ 


J 
d the amount of processing tin 


iD 


next state cf the system an 

the tauged job gets from the processor during the scjourn 
hime. We *hen repeat the procedure until ths accumuiacted 
rrocessing time for the tagged job meets the requirement of 


work required to comrletion. 


Bo. WCRK TIME 


It turns out to te especially convenient to measure +imes 
in terms of the amount of actual work or processing that has 
keen acccmplished on the tagged job. Let C(w) denote the 
number of jcbs unde@ragqeing service at a moment when exactly w 
units of processing kave been accomplished on the tagged 
Fob". We will also assume r(j) = 1 for all jj. Mve wakes OF 
accretion cf clock or response tim? at work time wis C(w): 
1f C(w) = 1 then the tagged job is alone and response 
(clock) time and work time advance at the same rate, while 
if C(w) = 17 the tagced job is accompanied by 16 others and 
17 units cf response time accrue for every single work tine 
unit. It fcllows that the response time for the tagged job 


requiring T units of processing time is simply 
T 
R(T) = fcqw dw. 
0 


The process Sccw) } is a birth and death process related to 
X{w). It has arrival and departure rates x: = Aj(N-j) and 
nV = Hj. All the sinulaticns described in this thesis will 
be dcne in work time. 

To oktain the ccnditional expected response time of a 
tegged job that requires T units of work time as in previous 
chapter, we remove tke conditicn that the tagged job entsred 
to find j-1 others initially present in the system by 
applying the same steady-state distribution cf the number of 
jobs in system, gq , found in (2.4). 
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Coeoetueo lt tTON FOR A 2Z-TERBEINAL SYSTEM 


1. Algorithm 


cnly twe terminals. Fach terminal submits jobs to be 
processed by the processor at rate A and the processor has 
a service rate of ub for jobs already present in the system. 
The service effort is allocated equally to all jobs present 
in the system at any time. MicsenOrS, if a job requiring T 
units of werk time enters te find the system empty, the 
arrival rate of the cther job will be A while there will be 
no departure. Similarly, if the job enters to find the other 
one already present in the system, there will be no arrival 
while the service rate will be Mr(2) 72. 

To transform the rates into the terms of work time 
we multirly them by 17r(1) and 2/r(2) respectively. Thus, 
the arrival rate, if any, becomes A/r(1) and the service 
rate, if any, is simply p. 

Fased on the above transition crates an algorithm to 
rerfcrm sSinulation for the conditional response time and 
eventually «he mean response time of a job requiring T units 
cf precessing time will be given as follows. 


a 
¢ Simulate response 


Algcrithm t¢ simulate =ume, Of a) job 
requiring ZT units of Erocessing time in a 2-terminal systen. 
Tet Ww, = amount of work time remains to accomplish 


for the tagced job. 
Co = amount of clock time accumulated towards 


the respense time of the tagged job. 
Bese . See We = f and cg = 0 


step 2: If the tagced job enters to find the system empty, 


ctherwise gc to step 3, generate an exponential time with 
parameter A/rt(1). Call this t'. 
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| 


eer sew, sce the conditional response times 
R, = Cy t Wo /T(1)- 
SO P 

EieheiG CNG, Set ° 


= =~ > 9 
Wo = Wort’, Cc 


GOVLO. Steps. 


0 Cau / 5 { 4) : 


Beep 3 : If the iob enters to find another one already 
present in the system, 1.2. j= 2, generate an expcnential 
with parameter pd - Galteaenas t'. 

ae) ert ee AW aa, 


Rote Cot 2Wo/t (2) - 


sat the response time 


STOP 
Ee ae SN ece. |S 
Wo = Wot", Co = Cot 2t" /r (2) - 


GO TO Step 2. 
Repeat the above precedure until we obtain the 
conditional response times for both cases. 
To find the mean response time we use the long-run 
distribution of the rumber cf jobs in system just after the 
tagged jcb enterad, 14 Cy 


qj a oT A je = cM pr (3) ’ 


Bere 7 = 1,2 and g,tg, = 1 and TT; is the staticnary 
distribution of the continucus time Markov chain, pate ye 

The mean respense time of a job requiring t units of 
work time is then 


Peni tee + Re dae 


where R,,R, are the ccnditional mean response times, given 
the tagged job enters to find the system empty and one job 
already present respectively, generated by the algorithn. 
The variance, skewness and kurtosis of the response 
time may also tke obtained by deriving the usual central 
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Boments. These calculations will be detailed in h 


(OD 

J 
m 
2s 
eat 


secticn. 
2. Numerical Results 


The numerical results shown in Table II are cktained 


from the outputs of a computer program written in FORTRAN 
= === = —" ——< 
| TABLE If | 
| Simulaticn results for 2 terminalis | 
| N= 1, We 2 | 
| Tine Nean Sstd.Dev. Skewness Kurtosis | 
| 0.1000 0. 1339 0.0452 Onoes9 -1.4462 ! 
| (.0006) | 
| C.2000 0.2648 O20,:8:5i0 Gn 500 shea 
(.0014) 

| C.3000 0.4031 CeZ229 Oro 459 -1.3044 
| (. 00:27 ) 
| 5.0100 0.6688 Oa oars 0.6477 alee Ss | 
| (.0047) 1 
| 1.0000 1.3669 C2292 00 0.6017 rl. O29 | 
| (. 0098) | 

2-00C0 2.6864 0 5:05 O56: 90 -~0.4870 | 
| (0187) | 
3.0000 4.0128 Om6 516 0.5184 a2 oS ’ 
| (. O201 1) | 
{ U0 CO 6.6863 0.8224 0.2974 =U io2 | 
| (.0269) | 
ee 2S cas: i se 
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USsang =he above algorithm to simulates the conditional 
respenss time, given the initial system state when the 
tagged jcb arrives. The steady state distribution is then 

aie 


used to ¢vaiuetes the unconditional moments, and measures of 


skewness and kurtosis for various work time requirements of 
the jceb. The values -tetween parentheses below the means are 
their cerresponding standard errors which depend on the 
Number of simulation replications. These outouts are evalu- 
ated ktased on 500 replicaticns of the response time for each 
initial ccendition. The r(j) 's are all assumed to be unity. 

Note that the moments obtained from the simulation 
agree well with those obtained by solving directly the 
system cf differential equations for the moments of the 
response times. This fact provides a check for the simula- 
tion. The simulated response times show diminishing values 
of skewness and kurtoSis as the required work time becomes 
large, again suggesting that there may be an increasingly 
accurate normal aprroximaticn eC the response eS 
meee ripution. 


D. SIMULATION FOR AH N-TERMINAL SYSTEM 
1. Algeritha 


Now, consider a more general computer system with 
cne processor and N terminals. As before, eéch terminal has 
a submissicn rate A +o the processor andthe processor 
processes each job to completicn with rate tH. =the alloca- 
tion cf processing time always follows the method of 
tPmocesscr-Sharing™. 

The simulaticn for response time of a job requiring 
T units cf processing time will be done under the condition 
that when the job arrives there are j terminals active, i.e. 
j jotks, including the tagged one that just arrived, are 
bkeing served by the CFU. Here, j can be 1,2,...,N. 


Zo 





Onder Processcr-Sharina scheduling, if j jobs ars 
present then in a shert time interval of length h the tage 
Jeiemegcca nity), ) Ubacs Of work done. Thus, if W(t) is the 
amount cf werk done cn tagged job by the time it has been in 
the system for t units of clock time, and if the number of 

oa) 7a) 
As long as W(*) is less than the required amount of 


jobs in system during this time tis j, then W(t) = 


work T fcr the taggec job, we will have to accumulate the 
amount ef werk done computed according to the number of jobs 
in system at that time. The conditional response time will 
te the clock time t for which the accumulated completed werk 
time Sweat: is equal to T. 

° Again it is convenient to measure time in terms of 
work time. The work time process, C(w), described in the 
rrevicus section is a birth-death process with trates 
A(N-j) 5 = rj and _ pj seeped, iy ' 5 3 = N- The response time 
is simply R(T) = $C (w) daw, ict (oe= 1. 

To simulate the work time process, we generate two 
exponential times with parameters r; and Hy respectively. 
The minimum of the two will indicate which event, arrival or 
departure, cccurs first. If an arrival occurs first and the 
accumulated completed work time is still less than «the 
requirement, T, the number of jobs being processed by thea 
Ca, 1.¢€. system state, will be j+1l. Likewise, when a 
departure cccurs first the number of jobs for next 
computation will be j-1. 

The above okservations allow us to construct an 
algorithm tc perform a simulation for conditional response 
time and lead eventually to the estimation of statistics for 
the response time of a job requiring T units of processing 
time as follows. 


Algorithm te Simulate response time of 


requiring Tunits cf proces 


ystem, given when the t 


= 
he tagged job begins processing there 
sing precessed. 


s 
are (j-1) jebs also 


Tox 





let w, = amount of work that remains +0 accomriish 


for the tagsed job. 
Co = amount PeGscck  cime accumulated towards 


the respcense time of the tagged job. 


(hn 
ct 
iW 
lo 
ja 
i" 
oO 
ct 
= 


—oaleecidme = \). 


peep 2: im "jJ—l, mises the jobienters to fird the system 
empty, otherwise GO TC step 3, generat? an exponential time 
with parameter ),. eaeenirs t*. 
2) Gece, soe the cCesponse time 
R, = Cg + We/T (1) 
STOP 
Piet’ <Sw., set 
Wo = Wo - t! 
Co = Co + t'/r(1) 
j= 2 
GO TO step 3. 


seep 3: If j = 2,3,---,N-1, generate two exponential tines 
with parameters X; ard H;- Call them t* and t" respectively. 
ay eeeemini(< tt) Saar, 
Ry = Co + Woj/T (3) 

SLUP « 
Dy men (ec) SW, S 
m Lf Bin(t*,t*) = t', ser 


set the response time 


Wo = Wo - t! 
Co = Co + t!9/T (3) 
Pot 


GO TO step 2 or 3 or 4 according to i. 
D2 )reeeante',t") = t", set 

Wo = Wo 7 t" 

Co = Cg + t" F/T (5) 

B= ors ii 

GO TO step 2 or 3 or 4&4 according to j.- 


el 








Step 4:3 iii Nk Le SO . all terminals become active, 
generate an expcnential time with parameter My Conl this 
+ 1 
ayer et Wor set the response time 
Bag? = Can Wg N/T {N) 
STOP 
bE) te =" See, Se. 
Wo = Wo 7 Dak 
Co = Gon. Soe (N) 
a =) Vien 


GO Temata 3- 
Bach run of the algorithm for fixed initial j gives 
a realizaticn of the ccnditional response tim2 of the tagged 
job given there are (j-1) other jobs in the system when tha 


tagged jck arrives fcr processing. 
2. Moments of Response Tin 


The simulaticn based on the above algorithm provides 
eae batch cf conditional response times of a tagged jek for 
each initial system state. 

Suppose we simulate a batch of size K for each 
conditional response time at initial system state j, for all 
a: Se 

foe Sceetc ocueiereatezed cond tional response 
temierg2 Ven the Initial condition is j ; (that is, the tagged 
job arrives when (j-1) other jobs are being processed), for 
He ey Z2ecee,N and K = 1, 2,6 ac,Ke 


Mathematically, we can use the averages of fener 


i lee Over cdeheDaten . £0 COMpute . empirical first, 
second, third and fcurth conditional moments given initial 
condition j respectively. Te compute the unconditional 
empirical mcments for response time of atagged jot that 
requires T units of frocessSing time, we can multiply the 


j-th conditional empirical moment by the steady-state 
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rrobability qj given by (2.4) 
j = ie eee, N The empirical mean, 
fourth central moments then are 
expansicns, 1.6. 
N 
A = a 
EC R(T) ] = 2834; = R 
J ml 
AN ee as 
Meret R (I) | = &R2q.-(R) 2 
jai vod 
IX Ne E 
E{ (R(T)-R)3] = & R3q. -3R S R2q. +2(R) 3 
eve oy ~ ot 
A N eae 
E[ (R(T)-R)*] = BRSq. -4R Z R3q. +6 (R) 2 
m= © se 
3 Jet 
= : | 
where R; = (= RY) / Ke foment = 1.253864 


However, in practice 


numerically unstable since the averages of the seccnd, 


and fourth moments over the batch may 


Hence, when we add or subtract these 


central mcments, Jt is sposcsible 


Froduces rcund-off errors which 


Therefore, we would rather rewrite 


terms cf cconditional expectations, 


Ef (R-E(R))" ] = ECEC (B-E(R)'1X(0) J], 


where X(0) is the number of jobs 


the CFU, 
arrives. 


including the tagged job, 


manner. 
One 


more numerically statle 
the 


kurtosis of the respcnse time. 


shown fcr computations 


ae 


the above 


the central 


This allows us to obtain the central moments 
Further 


variance, 


and Sum Over all 
Variances, coeerd . and 
computed by CoO wer 
(3a. 1) 
N 315 
Req. ~3(R) * 
24 (R)*, 
procedure can bs 
esata. 


be very large numbers. 
numbers to compute the 


ation 


kne 
be 


what comput 


may substantial. 


moments in 


1. 


2=1,2,3,4. (3:2) 


requesting processing by 


when the tagged job 
ine 
details will be 


skewness and 





Ke N 
let Rk, = (Z By) /K eric = eg + Een) the 


estimated ~ 


JN = 
Teepe | = 1228 ee). 
—" kas jot 
= 1S Seq- Ry +R, - 8) 2q; 
kot yol 
K 
= ao 2 . 7 ap = Re = 2 « «6 
SSE (Ry Ry) 2+ 2(Ry -Ra) (Ry -R) + (RAR) 2 Ja, 
K=1.3>1 
K 
K 
Since ERR i = 0, the estimated variances becomes 


N —? = 
werk) = 12 Zea att & (BR -B) 2aj- 


b=) ja 


This is the sampling version of the formula 
Var{R] = E{LVar (IX (0)) Jt+VarL FE (R{X(9)) J, 


which is a known general result that applies *o any random 
variable, R, that also depends upon another random variable, 
namely xX(0). The first compenent represents the overall 
variability of R for a fixed value of X(0), and the second 
component represents the variability of R due to the 


variakility in xX(0). 


To compute the estimated skewness we first computes 
the third estimated central moment of the response time by 
rewriting it ina cenditional expectation form as we have 
done for the computation of the variance. The third central 
foment is derived = shown below. 


ep (Fi, -R) 3] = iF Sey -R)3q 


x Ket des 
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ve - 

BCR, 5) 9 J 

N a —_— eZ _— 

>f (Ry -R;) 343 (Rj, -R; ) 2 (Rs -R) 
\ 


+3 (Ray. -R5) (R; -R) 24 (BR; -P) 3 Bais 


We can rns timo eracCoOmponent On the right hand side, 


nce (3 Po, 7 = 5 Guntic 


Bt ( Bees) = 1Z 3 -R, ead 85 S03 oar )2(R, -2) gq: 


kas g?! K kev ott 


Sz, “B) 3q, 


jaz 


The méastre of skewness of respons? time is then 


/- =~ 
Skewness(Riy J] = Beret) Si 
(Var( Ry }) ¥ 


fe 
[rn 
}+ 
=| 
jc 
to 
eu 

t 
im 
{Qu 
eas 
{Gc 
ce 
Oo 
wn 
lr 
{wn 


We start by rewriting the expression for the 


estimated fourth central moment as follows: 


K\- 
x MF 


we 
BL (Ry -R) 4) (Ey -R) #4; 


Mz - 


x 
ra 
r 

Zl 

ee 


zMx 


g7F% 4S 


*\- 


G. 
a 
cod 


Mx 
& Mz 


(Ray Ry) 444 (Ry - Ry) 3 (Rj -R) 


Al- 
x 


+6 (Ry, -R, ) 2 (R, -R) 244 (Ri, -B,) (By -R)3 
+ (RB; - R)¢ Iq) - 
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We then simplify, as before, by removing the 4-th term on 
the right hand side. Hence, 


a = K WN es: 
E{ (Ry -R) *] = 1% De eee ele HS DR, -R,)3(R, -B) q, 
x kev 421 K kava 


65 308 Sle 7B) 2 CR -R) 2a 2 -2) 4q, 
x ket = 


Therefore, the measure of kurtosis is 


A , 7 Sas 
Kurtoss[ Ry | Sete ye) Cs 
(Var( Ry }) 


Ssnce for the normal distribution zhe kurtosis has the value 
3, we then substract. = from the kurtosis computed above. Sc 
the new value of kurtosis will be 0 when the distributicn of 


response time has the normal degree of kurtosis. 


6. Standard Error of the Mean Response Time 


In crder to assess the accuracy of the simulated 
average response time, we may comput? a standard erroz for 
the mean response time of a tagged job that requires T units 
cf werk time from our batches of Simulated conditional 
response times. 


We have derived previcusly that the mean response 


time is 
~ 
/™ ae 
ELR;,] = R= 2B ASR a; 
uk jet 15, 


Now, we apply a property of the variance function by consid- 
ering q and K fixed constants and noting that the Rik are 
obtained frem independent realizations. This allows us to 
write the estimated variance of the expected response time 
as fellows : 
N 

an. 
Aj ass K 


™~ A» IK aia 
VarCEC Ry J) = 2(Ri Bi) 452 
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To ottain the standard error cf the mean respecnse time we 


A» ~ 
Simply take the square root of VarC ECR yy J ]- 


7, Numerical Results 


Tables III, IV and V show the outputs from @ simula- 
tion program based cn the previously described algorithm. 


The number of replications for each initial condition is 


| y 
TABLE III 

| Numerical results from simulation | 
| N = 10, AX = 15, PW = 100 | 
» Zane Mean Variance Skewness Kurtosis | 
| C.0100 -040368 .000295 .119865 -.715162 : 
| (. 000177) | 
| 0.0250 -100510 .001405 .088115 <-.712345 | 
| (. 000503) | 
| 0.0375 sisters .00@508 .0715783 -.587745 | 
| (. 0007 33) | 
0.C500 -201580 .003865  -.046461 -.572722 
| (. 000952) | 
0.0625 ~251340 .004953 -.012371 ~.510326 
| (. 001092) | 
0.10C0 ~402916 .008842 -.092651 -.378387 | 

(. 001510) | 

| Cc Sco -607487 .013806 -.061429 -.186003 | 
| (.001935) | 
| 0.2000 ~806743 .019704 -.023100 -.130595 
| (. 002303) | 
a ee eis. pean 


a 


if 








WO | 
| TABLE IV 
| Numerical results from simulation | 
| N = 10, X= 25, = 100 = 
| Tims Mean Variance Skewness Kurtosis | 
: O270'100 060567 200025 SoS «fa 5 050 | 
| ne | 
| rec 20 215069 2 -000992 Sent GOS |e OS ot | 
| (. 000467) | 

W257 50 ~ 226364 2007 5 woo? HS. sae45 0 | 
| a2 | 

7.0500 - 30240 8 ~90 2466 =——oozs 00 «38645 | 
| (.000791) ! 
| Ore Wierz 5 - 378130 200 S766 —souoo30 ~.527 906 | 
| (20009 25) | 
| 0.7000 603584 700 560m seosez?? §. 509166 | 
(. 00171228) | 
EE EE 


S00. The numbers below the means are their standard errors. 
We can see that tke results in the tables indicate a 
somewhat symmetric distribution for the response time 
(skewnesses are very small), and the kurtosises de not 
strongly indicate nen-normality they decrease towards 
Zero), especially when the processing time requirement 
becomes large. The kurtosis values in table III seem to 
indicate smaller tails than these in tabie IV, and thse skew- 
ness valves indicate that we have a more symmetric distribku- 
tion in the case A = 15 than when XA = 25 for an equal 
M = 100. In the next chapter some normal approximations to 
the respcenss time distribution wili be described. 
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TABLE V 


Numerical results from Simulation 


Time Mean Variance Skewness Kurtosis 

0.0100 .063016 .001011 .346171 -.463719 
(. 000202) 

0.0250 158884 .005204 .213798 -.555656 
(. 000627) 

0.0375 2229094 .010175 .170748 -.541305 
(. 001008) 

0.0500 2316021 .015985 .159817 -.565249 
(. 001354) 

0.0625 239478C .021523 .117365 -.496104 
(.00 16 €€) 


_ aE Ob a eee SP ESE a SE ae SPD ae aR PP PP ae ee EP ae eee EE Ee Ee ee ee ee ee ee See eee eee eee ee eee ee eee eee ee ee ee oP 


N = 25, A = 10, = 100 


Tine Mean Variance Skewness Kurtosis 

0.2101 00 -149847 .000659 -.527812 .471482 
(. 000252) 

0.0250 2375932 .002529 -.644450 .428531 
(. 000588) 

0.0375 -562173 .004470 -.667446 .366822 
(. 0008 34) 

C.0560 -750556 .006467 -.641420 .314407 
(. 001020) 

0.0625 .938072 .008259 -.593361 .275313 
(.001165) 


a ae ale 4b ae OSS ee, ae ee: ee a ae eee ae i a emer es eee 


—-. ...] 
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IV. NORMAL AFPROXISATION F 
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A. INTRCDUCTION 


The Markov assumptions we make on the processor-sharing 
system alicw us to infer that the distribution of tha 
response time of a tagged job that requires T units of 
processing time may be approximated by the normal distribu- 
tion when Tis large and/for when the system is in heavv 
feat iic. 

Iwo metheds are used to argue the approximate normality 

a 


of the distribution cf response time. Ona is based cn th 


a 


Central Limit Theorem for additive functionals cf a kirt 
and death process, and the cther follows from a heavy- 
Meret eac dittustcn approximaticn of the birth and death 
process. The formulas tc compute the approximate mean and 
variance will be described. More details concerning the 
analytic fcrm of tke appreximations are given by Gaver, 
Jacobs and Latouche. (Ref. 3] The approximations will be 
compared with the results freer simulation to study the¢ir 
@eeuracy. 

AS mentioned -befcre, the Central Limit Theorem fer addi- 
tive functicnals of VWarkov processes allows the ccnclusion 
that the accumulated work acccmplished in time t' of a job 
requiring a large amcunt of processing time, T, is approxi- 
mately normally distributed. This in turn allows’ the 
conclusicn that the corresponding response time is also 
approximately normally distributed. Hence, we will start by 
considering W(t'), the total work expended by the computer 
on the tagged job by clock time t' after its arrival, given 
that the tagged job requires exactly T time-units of work 
for ccmprletion. 

It ig cbserved that if when a job arrives there are 
X(0) = j customers, including the new arrival, present in 


the system for processing, then: 


oi6 





t 
ic) = J (200) 7%, (0)) da met) = 3 2), (re) 
0 


where Ke (%) is the number of jobs at the service stag at 
Greck tire t. 

Maem ths an aipemepriate central limit theorem for W(t) 
Bae De established by using results for finite birth-and- 
death models [{Ref. @j or ky making use of the theory: of 
cumulative processes. [{Ref. 9] We note here that the latter 
develcrcment of the central limit theorem is adaptable to 
models mcre general than the simple birth-and-death process. 

In «he case in which the system is in oe traffic we 
can approximate 4% X,(t);t20 3 by an OE Veen aUnlendeck 
process, s¢e¢ e.g. Iglehart. [Ref. 10] The process W(t') is 
then approximated by an integral of an Ornstein-Uhléenbeck 
process. A normal approximation for the response time 
@estribution then follows. 

After deriving scme formulas for the above two methods 
cf aprroximation, we will make comparison for goodness of 
fit of those approximations to simulation data. 


E. AFPRCXIMATION BY CENTRAL LIMIT THEOREMS 
1. A Central Linit 


Under Markov assumpticns throughout the processcr- 
sharing system, the number cf jobs being at service stage at 
eleck time <, Macy peas a finwte ergodic stationary time 
reversible Markov chain. 

if we define a functicn £ (X(%)) = xr(X (t))/X(t), then 
as outlined in Keilscn [Ref. 8] the process in equation 
(4.1) may fre proven ina variety of ways to be asymptoti- 
Coieeyvoenerm@at in distripution for large t. Hence, the 
accumulated accomplished work-time at time ee w(t) , 
Satisfies a central limit theorem. 
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Mmccccmmmomappevenrhe central limit theorem derived 
ty Keilscn we will have to redefine the infinitesimal gener- 
ator according to the number of jobs in system withcut 
Breluding the tagged cne. 


The relevant generator is now 


X'(t) = fj ---> X! (tth) 


Seale, or} (4.2) 


---> X! (+h) 


iF 
de 

( 
oad 


> Mi bi/ (j+1) +0 (h) 


-~--> X' (t+th) 


H 
Ld. 
ee 


= Nae M437 ata to: (bh), 


mom gy = 0, 1,-0.2, N'=N-1. 
The process in equation (4.1) then becomes 
t 


W(t) = Ve(xecuy) ay, (4.3) 
O 


meme ©£(kX' (uj) = c(X*' (uj) +1) / (X! (u)t+1). 
We can now express a statement for the central limit 


theorem ese: 


w(t) - zt ---> N(0,1), ast --->> © (4.4) 
oft 


foece the Constants j and (qe alte SUCh that 


4 


Ni 

5 = 26()T (4.5) 
J2O 

G2 = 2£(0)-£(1),2--,£(N") ] [T, zZf£(0)) (4.6) 
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Peeie2 Seong a Meertxy defined below 


2 = ({I-AtLJ™' - Lost, (7) 


ence ak identity Wace. x, and Lb is an N by N matrix whose 
° * ° ‘ = ° 
rows are steady-state probabilities, Lm, 4, OMe Pee ano 


death prteccess with rates in (4.2), 1.28. 
' , } 
Ee TT, e ¢ «6 TI Ay’ 


: ’ 
Tele 2. Thy (4.8) 


tt 
i 


creat enCTT A 


6 \ N 


and Ais a matrix defined as fcllows : 
Peel 


hd 
( 


5 ~ r,/78 ’ hoo = 1-dr>V% 


Asie = Ag 7ke Apgas = Pia(G7 4178 ARS = 1 ZS 4.9) 


a 
ul 


Mert Bien G/ (24), 8 = max 3, 


J 
(Yo = Aye Wy = My(k'7N)), (again with rates in (4.2)). 


2mit Theorem 


[Ih 


4 


or response time, R(T) 


The accumulated acccmplished work-time at clock time 
t' of atagged job requiring T units of processing time, 
Wet!) , increases in random straight-line segments from 
W(0) = O until W(t") = T. The value of t' hsre will ke the 
respense time of the job, R(T). It is the first-passage time 


to the required work time. So we can state the following: 


3 





Paictpee< fT) = P (RAT) o> t*) - Ee aon 


From equaeticns (4.4) and (4.10), we can derive, by the same 


standard argument of renewal theory as given in Karlin and 
meyebotefret. Way =pp. 208-209, the following statement of 42 
central limit theorem for response time, Ree): ; ae ous 


meprcaches infinity, i.e. 


pee a (275) 75> NW,1), as T ---> ©. . (Geeak) 


rene 


C. APPROXIMATION BY HEEAVY TRAFFIC ANALYSIS 


Ncw, we consider the processor-sharing system as 2 
Machine refair model in which A; = X(N-3)- Mi =H, and 
r(j) = 1 for all j's. Let N the number of terminals re large 
ema the traffic intensity Qs which is defined as the ratio 
of the expected service time tc N times the expected inter- 
arrival tine, be a fixed value less than one, 1.2. 
@ = B/ (XD - 

Under the above ccnditicns, when a job requiring T wnits 
cf precessing time enters the system it will ke processed in 
the company of many cthers. This indicates a@ system with 
heavy traffic situation. As mentioned earlier in the intro- 
ducticn section we may apply the properties of the limiting 
diffusicn process developed by Iglehart. [Ref. 10] 

Therefcre, 1£ X_(%) is the number of jobs at the 
Frocessing stage atclock time it, then X¢(t) can be 
approximated by a diffusion precess as follows: 


X(t) = Na(t) +{N .V(t), (4.12) 
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i jo cmdeascENoOn? Stic function of -ime and when * 
Peamenigmeo intense it has a2 finite 2imit, i.e. 2 (co) = 1-2, 
and $¥ (+t) Si, For the present model, A Olid me eNa.ae 
Ornstein-Uhliesnbdeck prccess. 
The accumulated ccmpleted work-time by fixed time =" of 
a tagged job requirin Taiipees Of processin time in 
equation (4.1) becomes 
. : 
W(t") = J au7x. cu) = Jauy (vaca) + fFex (a). (4.13) 
0 O 
Next, we apply the approximation and expand the expres- 


sion in (4.13) to seccnd order terms in N; then assume that 


the taaqged job arrives when the system is in steady state so 


that we can use the finite limit of a (0) =e ine Lace 16 = 


a(u). The arpproximaticn is now: 


t 


+" 
mace) = J any cvca- pn) : (AA/EN (1p) 12) | v(uy au, (4. 14) 
O O 


moO S t* < t. 


Hence, the expectation of accumulated amount of work- 


time ccmpleted on the taaged job is approximately 
t*/(N(1-Q)) ¢ and the actual distribution of total work dcne 
is alse approximately Gaussian (inteqral One an 


OCrnstéin-Uhlenbeck process), where the Gaussian prorerty 
results from the assumption of many accompanying jors, and 
not necessarily because the tagged job requires a leng 
rrocessirg time. 

As ** approaches infinity, we can evaluate (4.14) and 
show that a) Y(u)duj = 0, and vart '¥(u) au} = (2p/ (NAP) ) T- 
So the normal approximation to accumulated completed work 
time, W(t'), has the farameters: 


F = UN(H—) 6 G2 = ZH/( AZLN(1-—) 19) - 
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By réascning the sam? way as *0 obtain a central limit 
theorem for response time from a known approximated nermal 
distribution of accumulated completed work-time for a taqced 
job requiring T units of processing time, we can derive a 
normal appreximation for tke distribution of respense time 


as having the parameters as follow: 


BC R(1) 7 = N(1-p)T, (4. 15) 


Var(R(T)} = 2T/ ( APN (1-0) - (4.16) 


We can also improve the value cf the variance by using the 


formula telow: 


Var(R(1) ] = [2HT/ (Xs) Wee (Cone) 


| -ST 
where s = AN-B andf§f=1-[{(1-¢e. )/(s?)}. 


De. CCMPARISON TO SIFULATION DATA 


We will consider some particular simulation results from 
a system consisting cf 10 terminals and one processor in the 
cases of light traffic and heavy traffic to make comparison 
to nermal approximaticns described in previous section. 

First we can use the measures of skewness and kurtosis 
cf the response time resulting from simulation data _ to 
roughly assess the degree of normality or non-normality. We 
know that if a distribution is symmetric its skewness will 
tke zere and for the normal distribution its kurtosis has the 


value 3. Hcewever, eéevern if a distribution of the simulated 
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response time has tke measures of skewness and kuz*osis 


Bigcemtc tieee Of times nOmmal distribution, it does net imply 
comes chee dastributaicn is necessarily normal. ie Moya, 
suggests that a normal distribution may be a reasonatle 


mOproOximaticn. 

Secondly, es wouwne. compute the empirical distributicn cf 
the simulation data at various quantiles of the nermal 
distribution whose mean and variance are approximated by 
either central limit theorem or heavy traffic analysis 
(limiting diffusion). For example, to compare the one-tenth 
Gguantile cf the aprfroximated normal distributions to the 
Simulated response time at the initial systen state j, we 


first ccrpute 
Rio = Approx. Mean + (Approx. Std. Dev.) (Za) 


where Zio is the (1/710)th quantile of the standard nermal 
sas tribution. We then determine the conditional relative 
frequency of the simulated resfonse time, given the initial 


system state is j, as 


ref. 3 (Bo) = No. of simulated response time < Ruto. 





Nc. of replications 


we then use the distribution di in (2.4) to remove the 
initial state condition. Hence, we obtain the one-tenth 


guantile of the simulated response time (corresponded tc the 
normal appreximations). We may also compute the standard 
error of the estimated quantile by taking the square root of 
A N A A 
Var (EF) = (249,*p; (-P)) /K, 
jel d 
where E; = ref. 5 (Ro) and K = number of replications used in 


Simulaticn. 
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New consider a system of 10 terminals with arriva 
A= 15 and service rate HM = 100. The expected numb 


active terminals is, therefcre, 
FAL X (4) ] = 10(1- (1007150)) = 3.3333 , 


Mee n Indicates a light traffic situation. 


TABLE VI 
Ccmparison of Simulation Data to Normal Quantiles 


N =10,A = 15, U = 100 


HTA {I mOUU s0Z2d=eezoo «047 .934 .987 1.0 
an PVCs Foa250 62605 4893 .972 1.0 


| Time Normal .10 -25 -309 -25 -30 -93 -39 
| Cecnoo CLT 0.0 .123 .503 .887 .995 1.0 1.0 
| Ps 

| Hon (050 0.0 2261 .948 41.0 1.0 1.0 
| HTA II 0.0 .102 .361 .677 .909 .972 .999 
| 

im 0.0250 CLT ~042 .211 .487 .806 .967 .994 1.0 
| HTA I 0.0 .012 .334 .832 .996 1.0 1.0 
BTA II 0.0 .094 .334 .668 .914 .977 1.0 
| Se Os00 “err ~C&3 .232 .486 .766 .938 .977 .999 
| He rt  O.0 .032 .300 .737 <970 .998 1.0 
i Dewees 067s 300 2.649 «916 .973 1.0 
| C0750 CLT .082 .241 .494 .766 .922 .967 .997 
| ETA I 0.0 .033 .266 .688 .948 .991 1.0 
| HATA II .001 .049 .266 .630 .900 .969 1.0 
| Cealcoo Cur FCsvemeis -48s .761 -919° .969 .998 
| 


GO) 


+! 


(D 


Slee ee es a ee 


| 
| 
| 





—  @#«. — a. 
TABLE VII | 
| Ccemparison of Simulation Data to Normal Quantiles 
! 
N= 10, A = 25, WW = 100 
| | 
mame Normal .10 .25 -50 .J3 .20 .95 .39 
we0100 CLT moe os) 6 407° 4854 .997 1.0 lee | 
| MTA I mus eee .o54 .826 .990 1.0 le | 
| 
| Hike Tee@omae oo. .u3u . 721 .911 .975 1.0 | 
| O02 5SONCLT O77 .191 .453 .804 .975 .998 1.0 | 
ETA I ~084 .192 .438 .774 .955 .993 1.0 | 
| Prato toe e216 2.438 .729 .917 .976 1.0 
| O20500 “CLT wore eo oe 4od 6 1/69 4.951.990 1.0 | 
HTA I wee ec esa 72605927 1.979 1.0 | 
Pinte meen tCG e22055.4395 .707 .909 .965 .998 | 
| ORG 50 SCLT 093 .198 .447 .754 .936 .983 1.0 
ETA I -095 .193 .419 .713 .908 .967 .998 | 
| poet 5 102 200 419 701 .892 .959 .996 | 
| 
Oe OC Car rtUt ecco «4 fle. 756 .933 .976 .999 I 
| HTA I - 102 .226 .448 .711 .898 .958 .997 | 
Pia r . 106 «232 2.488 .708 .891 .953 .996 | 
ee hLhLhLCU tC 


Table III shows the values of skewness and nermal 
kurtcsis cf the response time of a tagged job for various 
work-time requirements. These values are quite small 
especially for large work-time requirements. It indicates 
the pessibility of a successful normal approximation. 

Table VI shows hew the simulated response time of a job 
requiring some processing time units is distriktuted 
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Sourerangqetc the Vaer0US quantiles ( .10, .25, 0, i, 
bY 
central limit theorem and the two heavy traffic (limiting 
diffusion) approaches (HTA I is (4.16), HTA II is (4.17)). 


We see that in this Situation the limitirg diffusion 


BOges9S, 299 ) Of the normal distributions approxinated 


eppreckematicn for normality of the distribution of resronse 
time is net so good an approach even for large amount of 
processing time requirement. This is easily anticinvated 
Since the system has only an ¢€xpected proportion Cit one. 
j@rrd of its 

theorem seems tec work pretty weil when the required work 


terminals active. However, the central limit 


time E¢eccmes very larce. 
Llet'sc 
and the service rate » = 100, with N = 10 ag 


onsider anew system with arrival cf 


expected number of active terminals is now 
EC X(t) ] = 10 (1-(100/250)) = 6, 


which indicates a moderately heavy traffic in the systen. 

Table IV provides us some gocd feelings about the meas- 
ures cf skewnesses and nermal kurtosises not being toc far 
from those cf actual normal distribution. Now, observe the 
table VII which shews the distribution of the simulated 
response time of a tagged job for various work time 
requirements comparing to the various quantiles of the 
hormal distributions approximated by central limit theorem 
and limiting diffusion approaches. 

We see that the limiting diffusion methods, in this 
case, wCcrk pretty well even for small work-time require- 
ments. The one with second formula (4.17) of approximating 
the variance works better than the other, because it 
provides a smaller standard deviation. The central limit 
theorem approximation still works very well when the werk 


time requirement becemes large. 
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Comparison of Simulation Data to Normal Quantiles 
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The above rasults for the two situa 
traffic and another moderately heavy traf 
confirmation that for a tagged job requiring large amount of 
processing time wee dis Ea bution ope resvonse time 
approaches, asymptotically, normality. In heavy «ratfic casa 
the liriting diffusicn approximations seem to work better 
than the central ilimit theorem approach. They do not werk 
well at allwhen the system has the terminals active léss 
men a half of its full capacity. 

Te ccenclude this chapter we show some more results in 
comparing the approximated normal quantiles to the values of 
Simulated resvonse times in weysie2 VIET for different 
(larger) number of terminals and various transitional rates. 
We can see easilv that as the number of termina 
the distribution of simulated response time approaches the 
approximated normals more rapidly. The heavy *raffic aprprex- 
imaticn seems tec work better, in thes2 cases, enan he 
central limit theorem approach, especially for small 
required work-time. This corresponds to our observation that 
the central limit theorem requires large work-time require- 
ment to tke a gocd aprroximaticn, while the only cendition 
that the heavy traffic approximations require is «+c havea 
certain amount of joks waiting tc be served at any instant. 
Finally, table IX shows the mean and standard deviation of 
the respcnse tim? comruted by the central limit thecrem and 
heavy traffic approximations for all the cases that we have 


keen studying the comparison of the quantiles. 
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| 
| TABLE IX | 
| Arprcox. Mean and Std. Dev. for One-Type Model | 
| | 
| N= 10, A = 15, M = 100 | 
CLT HTA I Ue 
| Tine Mean Std.Dev S2an Std-Dev. Mean sta-Dev. | 
| goto0 .0404 20322 7.0333 .0516 me s3 950236 | 
| 0.0250 .1010 .0510 .0833 .0816 .0833 .0535 | 
| Mose G 5.2019  .07z2 .1667 .1155 .1667 .0919 | 
0.0750 .3029 .0884 .2500 .1414 .2500 1216 | 
| Oeoue 24038 .1020 .3333 .1633 52333 1462 | 
nc | 
| N = 10, A = 25, BM = 100 | 
| 0.0100 .0605 .0245 0600 .0231 .0600 .0160 | 
! 0.0250 .1513 .0328& .1500 .0365 oo Os 14 
| @20500 .2027 .05428  .3000 .0516 -3000 .0481 | 
! 0.0750 .4540 .0672 .4#50C .0632 SES OOM OGIO | 
| 0.1000 .6053 .0776 .6000 .0730 -6000 .0706 | 
eee hh LULU | 
N = 25, A = 10, PB 100 | 
0.0100 .1500 .0390 .1500 .0365 SOO Skil 
Geos 0mece5>0 6.0617 33750 .0577 #£«+«.3750 .0497 | 
| 0.0375 . £625 Ggse) eco 258 2.0707 M025 wea | 
0.0500 .7500 .087Z .7500 .0816 -7500 .0760 | 
| 0.0625 .9375 .0975 .9375 .0913 .9375 .0863 | 
N = 40, A = 10, W = 100 | 
Gaomee . 3000  .0264 .3000 .0258 -3000 .0213 | 
| Mowe ees 0512 1. 1250 .0500 fen > 0 20077 | 
| @.0i625 1.8750 .0661 1.8750 .0646 1.8750 .0628 | 
| | 
\ Pee. il 





V. MODE 


Aw. INTRCDUCTION 


The processor-sharing model we have been describing so 
meeaed l= W=-th a Computer system having a single exponential 
service time distritrution only. To generalize it wea will 
consider a model of a computer system that consists of one 
Central Erecessing Unit, M terminals of type I and N termi- 
fees cf type II. Type I terminals submit jobs tc be 


Peeecssed by the precessoe at rat2 rA,, Which need an 


expected amcunt cf work 1/#P,. Likewise for type II termi- 
fey -L€ alcrivVal Late te the processor is X,, and «the 
expected work needed is 1/p,- Think times and amounts of 


work requested are tc te independent and exponential. 

This computer system may be viewed as the one having the 
ability to process two types of jobs. Each tyve of jcks must 
ke submitted from its correspending type of terminal. 

The expected resronse time of a tagged job that requires 
T units of processing time will be derived by the same 
approach we used for the one-type model. However, now, we 
have tc consider conditioning the given tagged job to be 
either type I or type II. The continuous time Markov-chain 
to be considered for the model is bivariate, with one vari- 
able being the number of jobs of types I being processed and 
the cther being tke number of jobs of type II tkeing 
fFrocessed. The response time of the tagged job will depend 
cn the numbers of both types of jobs in the system with it. 

Once we obtain the expected conditional response time 
given the initial condition of how many others as described 
above, we will have to remove those conditions by the 


steady-state distribution of the tagged job given it is of 
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ieee o> Cr type EI. The applacation of this steady stats 
dastributicn to the ccnditional response time wii 

to ccempute *he mean, standard deviation, and high 

ef response time of a tagged job that requires T units of 
work time, and hence the skewness and kurtosis as well. 
These values can help judge the goodness of the nermal 
aaprcXimaticn. 

Simulation will te used to generat? the conditional 
respense times as befere. The steady state distribution of 
the tagged job will be used te compute estimates of the 
mean, variance, skewness and kurtosis of response time for a 
tagged job requiring TI units cf processing time. Finally, we 
will describe an apfroach for normal approximations by the 
central limit theoren, and by a heavy traffic approach as 
for the one-type model, and study their accuracy through 


Srey iaticn. 


Be STEALY-STATE DISTRIBUTION 


When the system has two types of terminals to deal with, 
the direct derivation to find differential equations cr to 
apply the Laplace transforms for the conditional resfonse 
time of a tagged job is much more complicated than when the 
system has cnly one type of terminal. Since, we now have to 
consider net only tke conditicn that the job is tagged to 
find how many of them already present in the system Lut also 
the ccenditicn that tke job is of what type, I or II. 

However, under the processor-sharing concept with Markov 
assumpticns, exponential think times and work request times, 
it will not be tco hard to simulates the expected conditional 
respcense time, given that the job is one of the two types, 
and that it initially finds i-1 others of the same type and 
j of the other type present in the system. But then to 
remove these conditions we will have to apply the joint 
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steady-state distributions of the number of jobs present in 
the system just after the tagged job enters, given the 
tagged job is of specified type. The similar computaticn is 
performed for the case the taaqcad job is of anether +tycs. 
Let Mand N be the a of type I and type II térgri- 
nals respectively. If Re y(t ) is the number of type k jects at 
clock time = for k = on then, before a tagged jcb enters 
ene C(t Oru O,ls 2.0... gM and j = 0,1,2, «ee ,N, 
ke the system state i type I jobs and j type EL Jjots:ce2ng 
rrocesesed, the limiting distribution of the number of type I 


and type II jobs in the system is 


a @ e N— oN e e 
Wi,j) = dim Ef (xX, (t)-X,(¢)) = (i, 3) }.- 
t -» 00 
Consider the locel balance equations for type I tagged 
jobs: 

Net NZ 
TT(0,0) MD, = 4, TT1,0) 
Nw oS 
TT, 0) (M-1) Xd, = M 12,0) 
fw One 
tO CH =1yA, = Hy TT(i+1 60) (5.1) 


Viremeiiriy thas, tory = Of 12,...,M, the probability that 


the system is in state (i,0) is 


THis 0) = M(M-1)... (M-i+1) (7) TT(0,0). (S32) 
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Likewise, we can derive the same thing for type II «acged 
jobs: 
Tn Nd 
THO,0)NA, = BA THOe1) 
SNe NS 
TNO, 1) (N- NA, = BP, TO. 2) 
Tf _ “NZ ; 
Or 5) W- PA, = PATO. G+ (oimsy 
Mee tieeUen Thely Chac, ECE jf = O,1,.24,N, 
THC, 3) = N(N- 1)... (N-j+T) (A,7H,) TT(0,0) . (5.4) 
Other local balance equations are: 
LS 
TTH(i,0)N A, = 1 eee Mi. 1) 
b+ 
S~/S 
TH4.1) (N-1)A am ¢ 2) 
2: * he 
a 
TRG, 2) (N-2) - 2 wit. 3) (55) 
L+3 


ry 
THt1, 35) (HSaXx, Asp, Th, j+1). 


Lt +4 


Finally we can enya guess, from the equations (5.5), 
the steady-state distributicn of the total number of jcbs 
(of beth types) before the tagged job enters as, fOr 
Opes. em and j = 0,1,...,N, 


SS 





le : L+3 ; L 
TH. 5) = FS \ncr ne. cr-ien OX) 


J 


3 ~~ 
x ON (N=1) 226 (N-5+1) CAP T1000). (5-6) 


Mee diStra@cution alec satisfies the full balance equations. 


Giese stcady-state distribution of the number of jors in the 


system can be found by choosing TT(0 ,0) SC + Hae 
ma w Se 
= IT: 3) a Mere details can be found in Gever anda 
{sO JF 


Jacobs. [Ref 12] 

The steady state distribution of the antering tagged job 
mewmg Of tyre I and there being 1 jobs of type I and j jcbs 
Semcyps IL frocessing when the taqged job enters is 


ot, g7t) = k T{(1-1,)) (M-i)X., (oe) 


@na f£¢r type Il tagged job 


e ee. 
Mle, j,tt) = K 72, j-1) (N-S)A,- (5.3) 
MN 
where k is selected so that = < fecu@Meeie +g (1,.4,28n) ) = 1. 
26 jJ=90 


Peevey, scone ccnditaonal distribution of there being i 
jobs cf tyre I and j jobs of type II processing when the 
tagged jcb arrives given the tagged job is of type I is 


SN 
Or) = k ieee Tal) (ML) - al 


I 


MAN 
where k_ is chosen sc that Z21 Geaje= 1. 
1 lnO jzo 
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Oe ea CENTRAL LIBLT THEORES FOR THE RESPONSE TIME 


in thtewe2ction weewill present a central linit «theoren 


9) 


Bemeene GONd2.10Nnai distribution of the response times of 


emerged job requiring T unsts of work given i1t is 


sw 
= 
O 


particular type. Mere details can be found in Geve 
Meeobs. ([{Ref. 12] 

In what follows we will assume the tagged job is of 
meype Tt. Let K (?t) be the number of type I jobs (excluding 
the taaged one) being processed at work time t 3; that is 
when the tagged job has acquired «+ units of work. Then 
3 (x, (t),%,(*)); t20 4 PomaGOncrnuous time Markov chain with 


rates 
(i,j) ---> (i+1,5) : Cie (sey A (i+ 5+1) O<i<M-1 (5.10) 
(2,5) ---> (2,541) 2 (CHI A GF i+ 1), OS 5<N (sea) 
(i,j) ---> (i-1,3) : iy, See en) 
(i,35) ---> (i, j-1) : ip, 1S 4<N. (5.13) 


Similar arguments to those used in deriving (5.6) shew that 
ere litataing distribution TH(i, j) BOs the Markov’ “cGhein 
neweemommares) (>. 10) (5.11), (5.12) and (5.13) is of the forn 


THE. 5) =(“Jarn (H-2) 206 (HH E) OAH (5. 14) 
j 


XW (N-1)- 66 (N-J#1) CAPA)” TT0,0) 


3) 2) 





where 7T(0,0) is chosen so that SBM. 4) = 1. 
The response time of the tagged job which requires werk 
eo S 
+ 
R(T) = j [X,(u) +X, (up +ijdu, (Sa) 15) 


Oo 


memantegral of a funetion of a continuous *+ime Markov chain. 
Pomc Scott remecentoals limit theorem of Keilson 
fRef. 8] applies to give that as T ---> 0, there are 


GCometants mand g such that 


roe | TL 


Olt 


has an asymptotic ncrmal distribution with mean zero and 
variance one. In this case 


M-1IN 


B= 2 tisje) T,4) - (5. 16) 


LEO y2O 


The variance g? is ccmputed as follows. 
Let ¢€ be the infinitesimal generator of (X (7) 6X4 (*)) 
with system states in the order, 1.8. COO), 


eo eea@7 


Pee CMI lip 2 oe e (M151), woe , (0,N), weererti= Te Ns. Then 
C is an Mx(N4#1) square matrix having the form 


0 0 ODO «aa 


De week 0 8082... 


2 2. Ry 0 @ es ® 
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where the ron-zéro elements are the rates of those systen 
states. 

Let Y be the maximum absolute value of a diszgqonal 
emement ci ©, theme ry Uniformization of the chain we cktain 
a matrix A =I + eo 

eet ee=etat iesin( it) = TW(i,j), fe) = ew) and 
(1Tf); = = 11; (4) £5 (4)- Mecimeaceconding CO Keilson {Ref. 8] 
the contre Limes theorem variance tern, = en 
integral Sf, (a), X, (uy) du is as fcllows. 

° 


G2 = 2 mf, . Ore), «+--+. (MA 1 ZI £, (5. 17) 
t, 
Tn 
meee 2 = 1 (i-AtL] 9} - L, and 
at, &, . % 
L= |Mo Tee + Ty (5. 18) 
Th, TT, o @ @ WR 
If we define X= (Xo 6X, peee eX) = [aaron ole, then, 


suce [aA a= =1,0, we have 
fae OT yxX = f . (5. 19) 
g 
Polcweiveng both sides of (5.19) by Vf, and since TQ = 0 
and LX = TIX, we then have TTX = Wf. Thus (5.19) becomes 


“4 ex =f- (Tf), (52120) 


ez, be€teer, 
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OX = - ¥(f- Tf). [bs21) 


We Car then Solve equation (5.21) for K and thus the central 


limit theorem variance g* from equation (5.17). 


De HEAVY TRAPFIC APPROXIMATICN FOR THE RESPONSE TIME 


As before, consider a ccmputer system consisting of M 


eet =GrEMinale with arrival rate A,and service rate w,, 


and N type II terminals with arrival rate », and service 
rate MW, respectively. In this section we will present a 
heavy traffic approximation for the response time. More 


details of the appreximaticn are described by Gaver and 
Jacobs. [{Ref. 12] 

Let 1, =A,M, re ec toc — 7a asTo simplify the 
notaticns used in deriving heavy traffic mean and variance, 
we define the following expressions. 


FaeesteweesOnve LOD mM, which is the positive scluticn of 


ft 
the quadratic equaticn; 


OF “Lc +(PoL pt, ae) 
eee fio a? Pes Ps) *C%. 7 Ps | 
+m, *{L, STE a ob ar Je 


wien set ns = Cm, 7 (Ly (1-m,)) J = alates 


Let a, = -(L, (1-B,) -P,-L, (a, +m,) J 
a, = -(1L, (1-7,) ] 
b, = -{L,(c-mg) } 
Bb, = -(CLa(c-1,) -Pa-L2(m, +m) J 
q*= L, (1-m,) (m, +m) +p, 


ac 





oF = eae) (Tm) +B BS 


Now put SS, = (~ (ba +a, ) + f(b, +a) F-F (a, bl -ay by) 72 
, = (- (£, +a, )- f(by +a, ) 2-4 (a, B= a, by) ) /2 


(Note thet S andS are solutions to a@ quadratic equation). 


and S 


Also fut Pi = Meee) 7 (5, So) > B2, = wy. (Sa = Oe) 
Boos eaeee)7(5; =5,) ; By2 = 82/(5)-Se) 
nen) =n) 5 Ry5 = 78,7 (5, 7Soe) 
eee (5, -5) +: K,, = —b,/(S,-So) 
S3) pene 7 (e eets bb) = (b, -b,)/(S,5,) 
Xo = (a, ~4a,)/(a, 53-4,b,) = (8, 7-8,)/(S—95,) 
Pa, = (83,5, #17 (Sy-S,) § Bsa = (G5, #17 (SQ-S, ) 
Ra, = ~(%ySgtt)/(So-S,) 3 Kg2 = -(G2SQt1) /(Sp9-S,) 


Ncw, let X;(t) =numrezr of type i jobs, eee Ne 2 
Frocessing at work time t (tagg2d job not included). Suppcss 


the tagced job is of type I. Put 


VY, (t) = (Xi(4)-Mmp)7 OM. (5.23) 


THen, Ge) satisfies the following stochastic differential 


equations 


d¥, (t) = -a,¥, (t)+(-a,)¥, 4) +O, dw, (+), (3. 24) 


dy, (*) a (E) - (“by ) YT, (t) + Gdw, (t) ’ (Be 251) 


where W(t) and W(t) are independent Brewnian motions. 


a9 





The response time for the tagged job requiring = units 


cf precessing time is 


a 
Bett ) = \LX, (a) +e, (u) Jay (225 ) 
0 
+. 
= (cYF yy (wy oun, + 4% Yo (uj) +Mm, Jdu 
O 


= [Mm (m, +m, ) ee iF Jt Y, (0) +Y, (a) idl tis 


e 
Ieee 2 (¢) = SCY, (a +¥,(u) Jau. tinea 62) ( =) Satisfies the 


stochastic differential equaticn 
ae) = Y(t) +Y, (t) . (S27) 
The keavy traffic approximation for the distribution of 
R(t) is that R(t) has a normal distribution with mean 


(M(m,+m,) Jt 
and variance equal tc M(Var({Z(t) }). 
Pemmoeci@t ve dil expression for Varf[Z(t)] we will aprly 
results found ir Arneld [ Ref. 13] to solve the fcllcwing 
system of stochastic differential equations: 


dv (t) = AV(t) + Bd (t) ’ (7.25) 


Aw 


where W(t) =(Y¥, (t)-Y¥g(t),2(t) 1, dW(t) = (dW, (£), dW, (t) } 


and wal “a, 0 one 
A = -b, -b 2) 0 S, 


= 
—_ 
© 
O 
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Ppemecucwlcemy (SGe2.4) On page 130,Arnold states the solution 


Cimegqwascacon (5.28) as follows. 


- Fy s,t Sot Se 
Y(t) = [p,e° +k, Eo WN, (0) +E p27 +k, S' (0) (5. 29) 


c 
$oct-u) S, (t-u) 
{ote aie = JaW, (u) 


t 
S(t -u) $,¢€-4) 
Sot p.¢ *K, JdW,(u), 
O 


7 Sot st Sat st 
ae”) = [B2\e + hoy Ss Te meio 2° +K, 2 JX, (9) (520) 


t 
Solt -u Sit - 
(ote, e eK SEW) oe (uy) 
O 


So (t-u) S,(t-u) 
+K, e 


t 
font P28 > Jaws (a), 
0 


Sat t 
Zt) = ( ier bye +R, a! 1X, (0) (5.31) 


C s,t 
ue 30+ Base” aa ) J¥, (0) 


te 
$o(t-) S,(e-) 
le 8 + p,,° +K,, ¢ JaW, (a) 


t 
$o(t-u> $,Lb-4) 
+f, %3,+ 32° aieere Jaw, (u) . 
0 
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Sot $,t Se Sit . 
Beer eye 54,5 2060 and C, = %,,4+ Bie +K,,¢' - It 
tiem fclicws that 
Wer{ Z(2) ] = ce Var{ ¥, (0) J+ CeVer(Y, (9) ] 
+ 2C,C, (Cov[Y, (0) Yo (9) }) (5722) 
- t-U) 
2 Solt-u $,ct-uy 
OF -2 
- (t-“) m 
Sete an Bso +Ky,e “eau. 
oO 
cd 
as a) = 2 ae oles 24 
where ; Var[Y, ( )} = g, C Bi, ¢ Ky, jedu 
O 
i SoU $,Uu 
Sfp natn 
O 
of} 
So ' 
vartY,(0)] = a@2\t pare” +key 2 J2du 
o 806 
SoU S$, 
+ qc C B22° tiKsee? jedu 
O 


oo 
Sot S,;U Sed S,u4 
= sz Jt Pué +k, @ J[paie +K2,2 Jdu 
00 


4 jau 


I Ba.” *Ks 2° 


2 
+ q C Bi2° +k, 


O 


Sy 





Weel tutememe integral forms for Var[Y, (0) ], ean O. | an a 
Covf ¥,(0),Y,5 (0) J Mimentrerenlet 2.652). fc “SmMalins Lew +O soive 
Memeo tion anteg=ats to Gbtain a formula for the variance of 


the resrense time by heavy traffic approximation. 


E. SIMULATION 


Mee Al Gieceathm 


Fractically, the simulation for resvonse time in 3 
System with M type I terminals and N type II terminals is 
perftcrmed by applying th2 same procedure as for the one-tynoes 
model. This means that, for each initial system state of & 
tagged jcb cf each type, we generate the exponential arrival 
and departure times te be able to determine the next system 
state after a sojourn time. The simulation is for the work 
time process. 

If the minimum of the two exponential times is 
greater than or equal to the amount of work-time require- 
ment, we simply determine the response time by converting it 
into real time term. If not, we determine the next systen 
state, the work-time remaining to be accomplished and the 
accumulated clock time, and generate new exponential times. 
we rereat this procedure until we obtain a response time for 
each initial system state of a tagged job of each type. Once 
we obtain those conditional response times, we can easily 
calculate the mean, moments, and neasures of skewness and 
kurtosis by applyine the steady-stats distribution derived 
previcusly. 

We now descrike in detail how to simulate 4 condi- 
tional respcense time, given an initial system state is (i,j) 
and the tagged job is of type I, for a computer system 
consisting cf M type I terminals having ), and Pas arrival 
and requested work rates, and N type If terminals having A, 
and Mr 2s arrival and requested work rates respectively ; (i 


does net include the tagged job). 


Gre 





Lil soe mew OrR=- Lame the arrzt val reta of the next 


tagged jceb cf sither type is 


eet 1) oem (Nase 1 (2+ 5+1) , (5. 33) 


and the service completion rate of a job in system is 


pe ae + J}, 1]- (S239) 


We, then, generate two exponential times with parameters A 
and CP. Call them t' and t" respectively. 


Lignand( c* ae") 284 


or then set the response tims 


R(isj,T) = cy te (itjel), 


where Cg is the accumulated clock time (initially = 0) and 


Si is the amount cf work-time remains to be completed 
(initially = fT work time reguirement of the jcb). 
ie pomt' 2") < Wo, set 
= ¥, min (t*,t™) 
ee Co) tan (f°, c)(1+5+ 1). 
If min(t',t") = ~*, which means the system state is changed 


by an arrival, we generate two more axponential times with 
parameters (M- (i+1)) >, (it+j+ 1) and (N-5) A, (2 +54 1) respéc- 
tively. If the latter quantum of time is less than the 
previous one, i.@. an arrival of type II job cccurs first, 
the system state changes from (1i+1,3) SOmmetert less |) 
Otherwise, it changes from (i+1,j) to (i+2,j). In the 
contrary, if min(t',t") = t", which neans the system state 
is changed by a departures, we also generate two exponential 


times fer service time with parameters +, and JR, 
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Soeeecumyety. §r -neslatter -{S the minimum of the two, i.e. 
Pee PatLetne Of) type IT job occurs first, the sys tem stats 
Meee then Change Erem (1,j) to (i,j-1). Otherwise, ant 
Ganges from (1,]) te (1-1, 4). 

Fased on new syst2m state and new work-tine require- 
ment, we repeat tke same procedure until we obtain the 
corresrcnding response time. 

Note that we will have in total M times N+1 response 
times fer type I tagged job and Mtl times N for type IT 
tagged jek. We can tten apply the steady-state distributicns 
dy (2-5) and Gay (4, 3) +o compute th? msan, variance, skewness 
and kurtosis of the response time of a tagged jcb that 
requires T wnits of processing time as for the one <<yps 


model. 


2. Computation cf Simulated Mean, Variance, Ske2wness 


a. Simulated Mean Response Time 


Let Ry(i,j,T) and RL (i/j,1T) be the k-th sinu- 
lated ccnditional response time, given the initial systen 
state (tagged jok included) is (i,j) and the tagged job is 
cf type I and tyre II respectively. 

If we perform K replications for simulaticn of 
response time at each initial system state, then the mean 


respense time, given that the tagged job is of type I, is 


M oN 
FeO een dat) ag (1,5), (Bo 35) 
zo UFO 1 
_ K 
where R(i,j,I) = (ZR (ie5,T)) JK. 
bay 
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tk. Standard Frror cf the Mean Response Tine 


Under an assumption that Ry 
Eee Jy 12) are independent, we can compute t 
the estimated mean response time, given the mwaggecdmjeD 25 Of 
type I, based on the value obtained from equation (5.35) as 
follows. 


ernie a 2{2 > 3) (5S. 36) 
° lj 


I< 
where Var{R(i,j,1)] = CZ (Ri, 3,T)-R(4,5,)) 2 YK. 
=I 


To oktain the standard error we simply take the 


square reet of the expression (5.36). 
c. Variance of Simulated Response Time 


We start Ey computing higher céntral moments of 
respcense time, given the tagged job is of type l, by 
applying the same method as for one-type model. The second 
central roment which is also the variance of response time 


may be ccemputed, based on simulation data, as follows. 


Var{R ] 


K MA 
ZZPU Fy EeF/T)-F) 2g (4,5) 1 (5. 37) 


1 
K ke} i=0 j=0 


K MN e 
= V2220 (8, (2,5, 0) -FG,5,0)) ? 
kK k=; L=0 j20 


. (R (4, 3, T)-8) 2 )q (3,3). 


a 
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ee Skewness cf response time 
Pp 


The skewness of response time can be ccmruted by 


eremyeveal formula, i1.é. 
where My is the =e ememal TOMent and @ iS the Standard 


deviaticn of the resrense time. Won mOr. the catewleas: Ca OF 
the variance, 


KMN 
= 1 R eivl) ~ R)3 Beg cc 
Ps 12,22! (Ry (4.3 qe (ie) (5. 38) 


20 \=O0 
7 Ly By Jel) Saeleive segs ).) 3 
- 3 (Ry (i, JeT) -R (4, 5,0) ) 2 (R (36557) -8) 


+ (Ri, jer) -R) 2 lap. 3) ° 


€. Kurtcsis cf response tine 


The kurtosis is defined as U,= (py a) - 
where My is the fourth central moment and 9g is the standard 


deviation of the resrcense time. As before; 


Hy = 122 DR (iede0) - Bea Ged (5.39) 


ito j=O 


and Ey simple pclynomial expansion; 
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K MN _ 
Hy * IS]ZC (Ry Ge 5-0) -RG GTI) 9 
= 244 


ee canoe 1) —R (2, ,5)) 3 (RG. i,D)-2) 


k 


meen: , J, 1)) 2 (R(i, 4, 1) -8) 2 


+ 6 i 


* (Rie JeT)-B) 9 a, (4. 5) 


3. Nume 


Tables X and XI show the values of the means, heir 
standard ezrors (numbers between parentheses below the 
h 


cr 
iw 


means), standard deviations, skewness and kurtosis of 


tH 


empirical response times for a system consisting of 5 type 


terminals and 5 type II terminals. We also attach the value 


in 


cf the mean and the standard deviation computed by «he 
central limit theorem and the limiting diffusion in heavy 
treat fic. Those valués are computed for two cases. In the 
first case, Table xX, the arrival and departure rates are 
dr, = 30, tA, = 100 for type I jobs and 2 = 20, Pe, = 50 for 
type II jobs. But in the second case, Table XI, those rates 
are in teverse order, i.e. , = 7a My = 5) (0)F > = 30 and 
2 = 100. 

we see that the results shown in tables kK and XI 
indicate that in both cases the distribution of the empir- 
ical response time of atagged job that requires a fixed 
amount cf processing time is about the same. This indicates 
that the distributicn of the response time is almost 


independent of the jck type of the tagged job. Th: measures 
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TABLE X 


Simulation results for a two-type model 
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Skewness 
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TABLE XI 
Simulation Results for two-type model 
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TABLE XII 


Approx. Ncrmal VS Empirical quantiles 


= SS SS ee Se 


Cr 
HTA 
CLT 
HTA 
Cra 
HTA 
Chu 
HTA 
SE 
HTA 


Git 
RTA 
CLT 
HTA 
CLT 
ETA 
CLT 
HTA 
CLT 
HTA 


Xr, = 
- 10 


- O49 
2 114 
ey 
ae 
CES 
of ein 
- 0&8 
0 111 
Pp 6 SS 
Pet ie, 


2 O47 
= 209 
2062 
- 103 
aoa 
Seto 
205 
20 
087 
« 115 


304 Az = 20, Py 


suai 
mee 
2194 
0228 
e204 
0224 
eee 
226 
ce 
se 22u2 


2 


tn 


ao 
2221 
eal o 
mee 
Pip sl 
0229 
-181 
eae 
aS 
eS 


2457 
e442 
e462 
2442 
0455 
2432 
2459 
2434 
2465 
© 432 


© 433 
2443 
2425 
- 438 
© 429 
2447 
2417 
~ 431 
- 436 
0451 


71 


2869 
ic 
- 809 
ae. 
So 7 
oe 
ey 
696 
7S 
« /05 


- 838 
i 
e100 
OS 
a 4 
13 
2/24 
sOo3 
2745 
- i 


= 100, Loe 50 


228 
2919 
-978 
2916 
Oe 
5 Shee, 
I) 
oO 
945 
.889 


998 
0915 
IO 
5 ea 
ani 
Ne 
239 
ood 
SZ 
-898 


ao) 
oy 
oy 
aoe 
ao. 
neo. 
oo 
0 
- 986 
2 954 


oo 
ai9a1 
eet 
90 
oO 
oe) 
oe0 
—EGo 


lero) 
1.0 
tie 0 
oie 
eee 
eo 
1.0 
Soo 


120 
oe 


1.0 
moo 7 


| 


EE ee ES ei ED i oe ee ates in Se  —)  s 


unlieD que el emai ie eee RE eee Os Be, ee ee ei i eee ED eee SE mee ennai Ee SEED cee sei, nee en yee as ce eee at 


——— 


EES <2 SS LS SE EE EE ee a i a ns eee 





of skewness for these particular systems indicate that the 
feeb idutscn Of the empirical Gesponse time is a Ilitti¢ 
skewed to the left, andthe measures of kurtosis show 
thinner tails than it should be for the normal distribution 
even for a moderately large work-time requirement. We will, 
therefcre, observe the difference batween the normal quan- 
tiles, avproximated by the central limit theorem and by 
heavy traffic analysis, and the empirical response time to 
wake further judgement on the apvropriateness of the 
appreximaticns. 

We now compare the empirical response time at egach 
initial system state to the various quantiles of the ncermal 
distributions whose mean and variance are approximated by 
the central limit tkeorem and the heavy traffic analysis. 
The ccencérpt of the cemputations of the normal guantiles and 
the relative freguencies of the empirical respons time is 
the same asthat described for the one-type job model in 
chapter IV. The results shown in table XII indicate almost 
the same behavior in distribution of the empirical resrponse 
time fcr both cases. Again, this fact indicates that the 
distributicn of the response time is almcst independent of 
the jcb tyre of the tagged jon. Those values indicate that 
cur ncrmal approximations agree pretty well with the empir- 
ical response times cktained by Simulation, especially when 
the werk-time reguirement becomes large. They also agree 
with the results shown in tables X and XI where the measures 
cf skewness show a slightly left-skewed distribution and the 
measures of kurtosis (small negative values) indicate that 
the distribution of the empirical response times has a 
slightly flatter peak, fatter shoulders and thinner tails 
than the normal distribution. Table XII also shows that the 
empirical response time has a slightly higher mean than that 
appreximated by either central limit theorem or heavy 


traffic analysis. Hcwever, we can stili say that the 


V2 


' 





EE SE ED ED AS OS LS A A SS A SS TT A, AD AT aS LS A LD A A i OD a SS A A A A S ea Ah TS SAS AS RS RD OS iy 


CLT 
ATA 


O2OZS0 


eae 
ATA 


Over > 


CLT 
ATA 


C.0500 


ert 
ATA 


0.0625 


CLT 
BTA 


| 


Mean 


0.0724E 
(. 00014) 


O03 9 
0.07197 


0.18158 
(. 00035) 


0. 1809 € 
Of ess 


O2zrats 
{.00051) 


0.27144 
0. 26990 


0.86227 
(. 00064) 


Orcorg's 
Oe so 0c 


0.45220 
(. 00077) 


0.45241 
0.44985 


TABLE XIII 
Simulation results for a two-type nodel 


0.02491 
Oe 1326 


0.02902 


Oa aI 3C 
0.02708 


0.03886 


0.04823 
0.03590 


0.04723 


0.05569 


0.04326 


0.05566 


ORO e224 
0.04964 


73 


10, Py = 


Skewness 


O05 74375 


mie 101 O32 


Sera 


= eran 3 


-0.94014 


125, WW, = 25 


KUreeis2s 


=O oto Wy 


=-0. 76209 


=O Oo 


Ue O 275 


067 483i 


| 
| 
| 
| 
| 


= qe —ae ei ee ee eee OE ee aes oe eat co Codi ieee: income <200=8 sopgey -cd el) “ates eee) -ce cues: —) cee -e et ~ i e  es e d  teS  —e) me pt ee ee ee ms 


| 





| 
| 
| 
| 
| 
| 


ee ee ee ee ee 


HTA 


0.0250 


CLT 
HTA 


eon 2 


Cine 
HTA 


G.0500 


Cir 
ATA 


0.0625 


CLT 
HTA 


Results 


i 
= 
»” 


Sid. Dev. 


0.07169 
(.00014) 


O07 at 
Ong no 


Ong S77 
(. 00035) 


Oma a2 > 
Om 7 oe 


0.26904 
(.00049) 


0. 2688€E 
0. 26990 


0.35883 
(. 00062) 


0.35851 
Cao Jc 


0.44810 
(.00074) 


0.44814 
0.44983 


| 


TABLE XIV 
for a two-type model 


0.01335 


UeOrz 365 
Own € 


02050 2 


0.03746 
0.02708 


OeUs 75> 


O204USe7 
Oa G3 a0 


0.04561 


0.05297 


0.04326 


Deez 2 


Om 9.2 2 
0.04964 


VQ 


Kurtosis 


wo. 8 1156 


=O Go.) 


Oa 220519 


Slur 


Os 3) 2 


een 7 oa, 


soe 1G 


Oo 170 


=9.79135 


aa) GlOieZ 


eae! 


a ee ee 


a cer em cm i age les cl ee es I i ee eg Se ey Sc ee a i 


| 





oe 
| 
| 
| 


TABLE XV 


Approx. Nermal VS Empirical quantiles 


Time Normal -.10 .25 .50 .75 .90 .95 .99 
oo1to0 CLT SireueeeaG 2.845 .907 1.0 1.90 1.0 
HTA peice rene ts Ss £710 .905 .973 1.0 
o- 0250" .CcCLT Sesumeeetotemeus> «820 .993 ie Ob 6 tO 
HTA eee eee? 6.682 .901 .967 1.0 
we0S75 CLI ~-070 .18€ .428 .795 .981 .999 1.0 
HTA 2-113 .224 .414 .693 .899 .967 .999 
o.0500 “CLT Sooo se sas . 76 2.972 .998 1.0 
HTA © 110 .216 .422 .689 .896 .969 1.0 
o, Veiz> “CLT ~-0C@3 .200 .443 .765 .959 .997 1.0 
HTA ioe oe te GSS 6 B91 «2.958 1.0 


Time Normal .10 -25 -50 -75 -90 -.95 -99 
SONG ert NOgueedicesus7 2906 1.0 1.0 1.0 
ETA .110 .230 .458 .737 .931 .982 1.0 
0.€0250 CLT .061 .178 .451 .824 .991 .999 1.0 
BTA .121 .236 .464 .734 .930 .985 1.0 
(ees 55 (CLT .068 .191 .441 .790 .979 .999 1.0 
ETA mmweee 2s oe ets | 6 725 .927 .983 1.0 
pees G) Chet -081 .189 .430 .776 .975 .997 1.0 
ETA 2117 .232 .444 .729 .935 .985 1.0 
Cac625 > CLT muceme 19 e442 ~.772 .960 .995 1.0 
HTA 2121 .237 .456 .737 .923 .979 1.0 


om eee ee eo OE ee oe ee EM cote OD eee SE cone ee cee SE ee es ere EE ce ee ee ee OD ee ee, ee ee ee ee ee ce ee se ED ge ei cee pene i com ee ee a 


eas RR mag a ay MI aS Rf, eM I a) ee ee CY Ke OE gigi arent meme <sym 


| 
: 
L 


i 











Geeoeo buctten Cr the Empirical response ‘time in table XII is 
not far from being nermal especially when the work require- 


menz is large. The heavy traffic analysis gives a bstter 


iD 


apprcximaticn for small amcunt of work requirement, whil 
for large werk requirement the central limit theorem sShews a 
slightly better apprceximation. 

More results are shown in tables XIII and XIV for 
the case cf A, = 40, A, = 10, My = 125, A, = 25 and the cass 
of reverses rates respectively. The values in both tables do 
Meme differ much. This again indicates that the distribution 
ef the respcnse time is independent of the job «ype of the 


initial tagged job. The comparison of the empirical restvonse 


ct 


time tc the quantiles approximated by the central limi 
theorem and the limiting diffusicn in heavy traffic shown in 
table XV also suggests the same conclusion. Finally, tables 
XVI and XVII shew the results fora single job type model 
With the arrival and service rates from the average rates of 
the two job type model, i.e. N= 10, A = 25, W375. The 
values in those tables suggest that there mignt Fe scme 
particular cases where we can use the single job type model 
to aprroximate the two job type model. However, further 
investigaticn is not included in this thesis. | 
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A. SUMMARY 


Wem Orouges Up aecomecDp: Of proce2ssor-sharing and “hen 
described a model for processor-shared work-time allocation 
cf a computer system with one type of terminalis as a birth- 
death precess under Markov assumptions in the first chapter. 
In the second chapter, the mean response time of a tagged 
job that requires a fixed amount of processing time was 
derived kased on a ccnditicn that the tagged job anters to 
find an initial number of jobs present for processing. We 
also derived some higher moments of the response time by 
means cf differential equations which led +o 
computaticns of variance, skewness and kurtosis o 
distributicn. Scme numerical results obtained by solving 
those differential equations were then given as aneé 
for a system with two terminals and the arriva 
departure rates are 1and 2 respectively. 

A Simulation procedure for the response time, given an 
initial system state when the job was tagged to the systen, 
was descriked in Chapter III. wWe detailed an algorithm fora 
particular two-terminal system to check with the numerical 
method and generalized it for a computer system with a fixed 
number of terminals. Eased on those empirical response *+imes 
we showed how to compute the mean response time using the 
steady-state distribution. The standard errors of the mean 
response tire, the empirical moments and the measures of 
skewness and kurtosis of the response time were also derived 
in that chapter. In Chapter IV, we derived the formulas to 
compute the approximate normal mean and variance of the 


response time by ‘tke central limit theorem for additiva 


Teo 





wD 


MiifersOndrs OL a karth-déath process and also by th 


ime ing darrtusion aprreximat.on in heavy traffic situation. 
Eased cn those apfroximated means and variances, we 
explained hew to compute the various quantiles and compared 
them to the empirical response times obtained by simulaticn. 
The difference between the relative frequencies of the 


|- 


empirical data and the CDF of the approximated ncrma 
distributions indicates how well the normal distributicns 
agree tith the empirical respens¢e times. 

In Chapter V, we described a model for a computer system 
with two tyres of terminalis, again under Markov assumpticns 
for a Eivariate birth and death process. Only now we have to 
consider a condition that the tagged job is cone of the two 
types. The steady-state distribution was thén derived for 
this two-tyrpe system. This allows us to remove the condition 
of initial system state and leads to the computaticns of the 
wean response time cf atagged job that requires a fixed 
amount of processing time as fer one-type model. The cemrfu- 
taticns cf higher moments, thus the measures of skewness and 
kurtosis, were also explained as well as the standard error 
cf the méan. Twe approximations for a normal distributicn, 
one LEasé¢ on the central limit theorem fer additive func- 
tionals and the cther one based on the limiting diffusion in 
heavy traffic situation, were derived for this two-type 
model in the same way as for the ona-type model.Next, a 
procedure +o simulate the response time of a tagged job 
given an initial system state was described, and then based 
cn thes? empirical data we make comparison to the quantiles 
of the approximating normal distributicns to study the 
behavior of the empirical response times towards normality. 

Finally, for appendixes, we attached two programs 
writter in FORTRAN te perform the simulation of the response 
times as described in Chapters III and V respectively. The 


program in appendix A simulates the response time cf a 
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tagged jceb based on its work-time requirement, the size o 
the system and the arrival and departure rates for a single 
job type model. It computes the statistic 2lements of the 
empirical response times such as mean, variance, skewness 
and kurtosis. It alsc computes the quantiles of the nernal 
distributicns whose mean and variance are approxinated by 
the central limit theorem and heavy traffic analysis. I+ 
then ccmputes the relative frequencies of the empirical 
tesponse times ccempared to the quantiles. The standard ¢rror 
cf each relative frequency is also computed. This will help 
us tc mak? a judgement on the number of replicaticns needed 
for simulation. The rfrogram in appendix B provides the same 
things fcr a two-type job model. We only have to input the 
méan and variance for the approximation by the céntral limit 
theorem from an APL frogram due to a large number of matrix 


operations. 


E. CCNCLOSION 


The agreement of the normal approximations with the 
Simulaticn is satisfactory for both one-typ2 and two-typs 
models, especially in the cases of large work-time require- 
ment andycr under heavy traffic situations. The central 
limit <cheorem apprceximaticn always works when the work 
reguirement is large enough ¢ven if we don't have a heavy 
traffic situation. The appreximation base on the heavy 
traffic analysis seems, however, to work better, given that 
the system is under a heavy traffic situation, than the 
central limit theorem approach for small work requirement. 
E@eeoieeael2G0> tkaftic Situation, the limiting diffusion 
apprceximaticn does net work well at all. 

Based cn the values of skewness and kurtosis of the 
empirical response time, we have tried to use the Edgeworth 
expansion tc improve the normal distribution of the response 
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time approximated by «he central limit theoren. The 
distribution computed by the Edgeworth approximation is 
Close to the CLT normal distribution, and the results dc not 
indicat¢ significant improvement. 

In the two-type mcdel, we observe that the response time 
distrikuticn appears to be almest independent of the tyre of 
the tagged job. Tt is hard to make a judgement on 
behavior of the empirical response time based on the obkser- 
vaticrse cf the mean, variance, skewness and kurtosis, eéven 
though those values indicate that the distribution is well 
apprcximated by the normal. The Simulation for the tesronse 
time cf the *two-type model involves a lot of computations. 
For example, in a system consisting of 5 type I and 5 
type II terminals, we have tc consider 30 initial system 
states in tctal. If we want to rerform &@ Simulation with 500 
teplicaticns, it will involve at least 15,000 computations. 
This fact indicates that the round-off error may be substan- 
tial. Hcwever, the relative frequencies of the empirical 
response time comparing to the approximate normal quantiles 
cbtained fer the particular cases in this thesis are ina 
rather gcod agreement with the theoretical suggesticn. This 
means that, under tke Markov assumptions, the distribution 
of the respcense time of a tagged job that requires a fixed 
(large) amount of precessing time is approximately ncermal. 

Threrghcut this thesis we have been assuming that the 
distributions of the errival and service times are exponen- 
tial. However, all the computations done for both models can 
as well te extended to the case of general distribution of 
se€rvice time. It might be of interest to study the models 
and derive a method tc simulate the response time, since in 
the real world we would always have to deal with a more and 
more sophisticated ccmputer system which is capable of 
dealing with many kinds of work. So further study might be 


needed fcr such systen. 
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APPENDIX A 
SIMULATION PROGRAM FOR ONE-TYPE MODEL 
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